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Summary 

—  •  In  this  paper  the  fully  nonlinear  theory  of  finite  deformations  of  an 
elastic  solid  is  used  to  study  the  elastostatic  field  near  the  tip  of  a  crack 
The  special  elastic  materials  considered  are  such  that  the  differential  equa¬ 
tions  governing  the  equilibrium  fields  may  lose  ellipticity  in  the  presence 
of  sufficiently  severe  strains. 

The  first  problem  considered  involves  finite  anti-plane  shear  (Mode  III) 
deformations  of  a  cracked  incompressible  solid.  The  analysis  is  based  on  a 
direct  asymptotic  method,  in  contrast  to  earlier  approaches  which  have  de¬ 
pended  on  hodograph  procedures. 

The  second  problem  treated  is  that  of  plane  strain  of  a  compressible 
solid  containing  a  crack  under  tensile  (Mode  I)  loading  conditions.  The 
material  is  characterized  by  the  so-called  Blatz-Ko  elastic  potential.  Again 
the  analysis  involves  only  direct  local  considerations. 

For  both  the  Mode  III  and  Mode  I  problems,  the  loss  of  equilibrium 
ellipticity  results  in  the  appearance  of  curves  ("elastostatic  shocks")  is¬ 
suing  from  the  crack-tip  across  which  displacement  gradients  and  stresses 
are  discontinuous.  / 

★ 

The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an  in¬ 
vestigation  supported  by  Contract  N00014-75-C-0196  with  the  Office  of  Naval 
Research. 


INTRODUCTION 


Among  a  number  of  recent  papers  devoted  to  the  study  of  the  struc¬ 
ture  of  finite  elastostatic  fields  near  the  tip  of  a  crack,1  the  investi¬ 
gations  summarized  in  [4],  [5]  and  [6]  are  concerned  in  particular  with 
elastic  materials  whose  corresponding  equilibrium  equations  are  enable 
of  losing  ellipticity  at  sufficiently  severe  strains.  All  three  of  these 
papers  treat  the  so-called  "small-scale  nonlinear  crack  problem"  asso- 
caited  with  the  finite  anti-plane  shear  (Mode  III)  of  an  infinite  slab 
containing  a  crack  of  finite  length  and  deformed  to  a  state  of  simple 
shear  at  infinity. 

The  principal  feature  of  the  results  reported  in  [4],  [5]  and  [6] 
is  the  appearance  of  two  curves,  issuing  from  each  crack-tip  and  termi¬ 
nating  in  the  interior  of  the  body,  across  which  the  displacement  gradi¬ 
ent  and  the  stresses  are  discontinuous.  Such  "elastostatic  shocks"  have 
been  discussed  in  general  terms  elsewhere  [7]. 

The  analysis  in  [4],  [5]  and  [6]  depends  critically  on  the  fact  that 
the  finite  anti-plane  shear  problem  is  governed  by  a  second  order  quasi- 
1  Inear  partial  differential  equation  and  can  therefore  be  successfully 
treated  with  the  help  of  the  hodograph  transformation.  For  the  more  im¬ 
portant  problem  of  the  plane  deformation  of  a  cracked  slab  by  tensile 
loading  at  infinity  (the  Mode  I  problem)  the  hodograph  transformation  is 
not  available  because  the  associated  quasi-1 inear  system  of  differ¬ 
ential  equations  Is  of  the  fourth  order. 

<1 - 

See  [1],  [2],  [3]  for  reviews  of  recent  work  in  this  area. 


In  the  present  paper,  we  first  show  that  the  qualitative  features 
of  the  results  in  [4]  and  [5]  for  the  Mode  III  crack  problem  can  be  ob¬ 
tained  by  a  direct  local  asymptotic  analysis  which  makes  no  use  of  the 
hodograph  transformation.  The  advantage  of  the  procedure  used  here  lies 
in  its  applicability  to  the  plane  strain  Mode  I  problem  which  constitutes 
our  main  objective.  We  determine  the  qualitative  structure  of  the  crack- 
tip  field  in  the  latter  problem  for  the  so-called  Blatz-Ko  strain  energy 
density  introduced  in  [8]  in  connection  with  experiments  on  a  highly  com¬ 
pressible  rubber-like  material.  It  is  known  (see  [9])  that  this  material 
is  capable  of  a  loss  of  equilibrium  ellipticity  at  severe  deformations. 

In  Section  1  we  quote  pertinent  results  from  the  nonlinear  equilib¬ 
rium  theory  of  homogeneous  and  isotropic  elastic  solids  and  we  introduce 
the  special  deformations  and  materials  appropriate  to  this  study.  Sec¬ 
tion  2  begins  with  the  formulation  of  the  Mode  III  crack  problem  dealt 
with  in  [4].  We  obtain  the  asymptotic  representation  near  the  crack-tip 
of  a  number  of  solutions  of  the  displacement  equation  of  equilibrium 
valid  on  overlapping  domains.  The  final  solution  is  then  generated  by 
a  consistent  matching  across  two  syrnnetrically  located  elastostatic 
shocks.  In  Section  3  the  tension  crack  problem  is  treated  in  an  analo¬ 
gous  manner  and  we  find  the  corresponding  asymptotic  solutions  to  lead¬ 
ing  order.  Section  4  is  devoted  to  higher  order  considerations  and  a 
discussion  of  the  results. 
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1 .  PRELIMINARIES  FROM  NONLINEAR  ELASTOSTATICS 

In  this  section  we  present  a  brief  sumnary  of  the  equilibrium  the¬ 
ory  of  finitely  deformed,  homogeneous  and  isotropic  elastic  solids.  The 
two  special  deformations  and  materials  relevant  to  this  study  are  then 
discussed. 

Let  ft  be  ar.  open  region  occupied  by  the  interior  of  a  body  in  its 
undeformed  configuration  and  denote  by  £  the  position  vector  of  a  ma¬ 
terial  point  in  ft.  A  deformation  of  the  body,  indicated  by 

y  =  y(x)  =  x  +  u(x)  for  all  x€ft  ,  (1 .1 ) 1 

rv  /v  /v  /v  rw 

is  a  mappi,.;  of  ft  onto  a  domain  R*  in  which  £(x)  is  the  displace¬ 
ment  field.  We  assume  the  transformation  (1.1)  to  be  invertible. 

Let  £  be  the  deformation-gradient  tensor  field  associated  with 
the  mapping  (1.1)  and  J  its  Jacobian  determinant,  so  that 

F=vy  ,  J  =  detF>0  on  R  .  (1.2) 

'  For  an  incompressible  material  the  deformation  must  be  locally  volume 
preserving,  whence  J  =  1  on  ft.  Define  C  and  G  by 

C  =  FTF  ,  G=FFT  .  (1.3) 

r**  A* 

These  deformation  tensors  have  common  fundamental  scalar  invariants 
given  by 

Vetters  underlined  by  a  tilde  represent  three-dimensional  vectors  and 
tensors . 
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Il  =  tr(C)  ,  I2  =  ^[(trC)2- tr(C2)]  ,  I3  =  detC  .  (1.4) 

To  continue,  we  let  be  the  actual  (Cauchy)  stress  tensor  field 
on  ft*  and  a  the  corresponding  nominal  (Piola)  stress  field  on  ft; 
t  and  a  are  related  by 

t  =  J-1a  FT  ,  a  =  Jt(FT)_1  .  (1.5) 

/V  rw  iv  /v  rw  '  ' 

The  balance  of  linear  momentum,  in  the  absence  of  body  forces,  leads  to 
the  following  equivalent  alternative  forms  of  the  local  equilibrium 
equations: 


div  t  =  0  on  ft*  or  diva  =  0  on  ft  .  (1.6) 

Suppose  S  is  a  surface  in  ft  and  8*  its  deformation  image  in 
ft*.  Let  n  and  n*  be  the  respective  unit  normals  to  8  and  8*  so 
that  the  associated  nominal  and  true  surface  tractions  are  given  by 

s  =  on  on  8  ,  t  =  t  n*  on  8*  .  (1.7) 

^  /V  «■>>  /V  ^ 

It  can  be  shown  that 

£=0  on  8*  If  and  only  if  ss0  on  8  ;  (1.8) 

this  proposition  Is  Important  because  it  allows  the  boundary  condition 
on  a  traction  free  surface  8*  to  be  specified  on  the  undeformed  sur¬ 
face  8. 

Let  W  be  the  stored  energy  per  unit  undeformed  volume  character¬ 
istic  of  a  given  elastic  material.  For  compressible  materials 

O 

W»  H(I-j  .Ig.Ij)  and  the  corresponding  constitutive  relation  can  be  given 


i 
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in  either  of  the  following  equivalent  forms 


t=  e,G+ z.,G  +  i-l  , 


-i 


a.=  J£.<si  +  e3£‘  > 


(1  *  9 ) 1 
(1.10) 


where 


P  _  _2__  /  3W_  +  ,  3W_\  .^2_3W_ 

1  ”  ^3  V3I1  ,3I2/  2'/I5m2  ’ 


v  =  2  /r  — 

l2  31, 


(1.11) 


For  incompressible  materials,  W  =  W( I-j , Ig)  and  the  corresponding  con¬ 
stitutive  law  takes  the  form 


=  2 


* 

3W 


* 

3W 


|t-G+  --(1,1  -  6)G 

3 1 1  ~  3l2  1~  ~  ~ 


-P  1  , 


(1.12) 


or,  equivalently. 


=  2 


* 

dW 


3W 


|=-F  +  fy-(I,l  -G)F 
3l<j  ~  3 1 2  l~  ~  ~ 


(1.13) 


Here  P  is  the  arbitrary  hydrostatic  pressure  required  to  accommodate 
the  constraint  J*/T^=l. 

We  now  turn  to  two  special  classes  of  finite  deformations:  anti¬ 
plane  shear  and  plane  strain.  For  these  the  region  ft  occupied  by  the 
undeformed  body  is  taken  to  be  cylindrical  and  a  fixed  cartesian 

^1  stands  for  the  idem  tensor. 


$ 
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coordinate  frame  is  chosen  so  that  the  x^-axis  is  parallel  to  the  gen¬ 
erators  of  ft.  Let  £  be  the  cross-section  of  ft  in  the  plane  x^ *  0. 

To  begin,  we  assume  that  the  material  occupying  ft  is  incompress¬ 
ible.  A  deformation  (1.1)  on  ft  is  an  anti -plane  shear  if  it  is  of  the 
form 

ya  =  xa  .  y3  =  x3  +  u(xi,x2)  .  (1-14)1 

As  shown  in  detail  in  [10],  the  deformation  (1.14)  can  in  general  only 

★ 

be  sustained  by  materials  whose  strain  energy  density  W-W(I^,I2)  is 
of  suitably  restricted  form.  Here  we  consider  a  class  of  materials 
shown  in  [10]  to  be  compatible  with  (1.14)  for  which  W  depends  only  on 

V 

W-Hd,)  ,  W/(I1)>0  for  1^3  ,  (1  .lS)^ 

W(3)  =  0  ;  (1.16) 

*.  * 

here  W  is  the  derivative  of  W  with  respect  to  1^ .  As  indicated  in 
[4],  substituting  (1.14)  into  the  constitutive  law  for  incompressible 
materials  (1.12),  (1.13)  permits  the  reduction  of  the  field  equations  to 

[W'UJu  ]  -U  on  »  ,  (1  - 1 7 ) 3 

I  » LX  j  CL 


^reek  subscripts  take  the  range  1,2  while  Latin  subscripts  assume  the 
values  1,2,3.  Repeated  subscripts  are  sunned. 

2 

Note  that  I i *  3  in  the  undeformed  state. 

3 

Subscripts  preceded  by  a  comma  Indicate  partial  differentiation  with 
respect  to  the  corresponding  material  cartesian  coordinate. 
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wlth 

I,  =  3+  |vu| 2  ,  |vu|2  =  u  u  .  (1.18) 

I  I  Cl 

The  components  of  Cauchy  stress  are  given  by 

•  T33-2“'(1l)l'?u|2  .  0.19) 

To6  =  °  • 

An  elementary  solution  of  (1.17),  corresponding  to  a  state  of  sim¬ 
ple  shear  with  amount  of  shear  y,  is  given  by 

u(x^,x2)  =  yx2  for  all  (x^,x2)€fi  .  (1.20) 

For  this  deformation 

t23  =  t32  =  =  2W(3  +  y^)y  .  t33  =  t(y)y  (0^y<od)  .  (1.21) 

We  refer  to  the  graph  of  t(y)  vs.  y  as  the  response  curve  in  simple 
shear  for  the  material  at  hand.  Because  of  the  inequality  (1.15),  the 
displacement  equation  of  equilibrium  (1.17)  can  be  shown  to  be  elliptic 
at  a  solution  u  and  at  a  point  (x-|,x2)  if  and  only  if 

t'(y)  >0  ,  y=  |vu(x1,x2)|  .  (1.22)1 

Thus, the  condition  of  ellipticity  is  satisfied  if  the  slope  of  the  re¬ 
sponse  curve  in  simple  shear  is  positive  at  an  amount  of  shear  equal  to 
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the  magnitude  of  the  local  displacement  gradient  vuCx^.Xg). 

The  incompressible  material  to  be  considered  here  was  introduced  in 
[4];  its  response  function  t(y)  in  simple  shear  is  given  by 


(my 


t 


(yH 


MY 


>1/2 


for  O^ysl 
for  1  *y<oo 


(1.23) 


where  m  is  the  shear  modulus  at  infinitesimal  deformations.  The  re¬ 
sponse  curve  (1.23)  is  described  in  Fig .  1 .  The  elastic  potential  W(I-j) 
associated  with  this  material  reduces  to 


«(!-,)  = 


JdT-3) 

+2y(l1  -3)1/4 


(3*1^4)  , 

(4  s  1-j  <oo)  . 


(1.24) 


We  observe  for  this  material,  the  differential  equation  (1.17)  is 
elliptic  at  a  solution  u  and  at  a  point  (x^^)  if  |vu|<l;  it  is 
hyperbolic  if  |vu|  >1.  In  this  paper  we  shall  study  weak  solutions  u 
of  (1.17)  which  are  continuous  and  have  piecewise  continuous  first  and 
second  partial  derivatives  on  &.  Clearly,  these  continuity  require¬ 
ments  allow  for  the  possibility  of  finite  jump  discontinuities  in  vu 
as  well  as  in  the  stresses  across  curves  in  &.  Equilibrium  requires 
that,  across  such  a  curve  C,  the  axial  component  of  Piola  traction 
given  by 


(1.25) 


must  be  continuous.  In  (1.25)  au/3n  is  the  normal  derivative  of  u 
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on  C.  We  call  such  a  curve  an  "equilibrium  shock11. 

Let  us  now  suppose  the  material  occupying  the  cylindrical  region 
ft  to  be  compressible.  A  plane  deformation  of  ft  parallel  to  the 
(x-|  .XgJ-plane  is  described  by 


*a  =  X»tUa<Xl-X2)  • 
y3  =  x3  ‘ 


(1.26)1 


For  a  plane  deformation,  the  fundamental  scalar  invariants  are  related 
by  I2  =  I]  +  I3  "  1 »  so  t^e  strain  energy  density  W(I^,I2,I3)  for 
compressible  materials  can  be  written  as 

W(IrI2,I3)-fi(I,J)  ,  (1.27) 


where 


1=^-1  ,  0  =  /Ij  . 


(1.28) 


Equations  (1 .26)-(l .28),  in  conjunction  with  (1.9)-(1.11)  and  (1.6), 
provide  the  coordinate  equations  of  equilibrium: 


2||v2ya+2 


JT 


9130  ,6) 


ra,8 


e8YGap  yp , 


=  0 


(1.29)2 


'One  may  consult  [12]  for  a  detailed  treatment  of  plane-strain  deforma¬ 
tions. 

2 

eag  are  the  components  of  the  two-dimensional  altenator. 
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together  with  the  in  plane  stress  components 


.2  3W„  w  .  3W  , 

Tae"Tea"T  3Tya,pyB,P  +  3J5aB  00  &  * 


o  o 

—  O  w  X  £> 

CTaB  3lya,B  3J  £By  apyp  ,y 


(1.30) 1 

(1.31) 


Furthermore, 

J  =  yl  ,ly2,2 -yl  ,2y2,1  ’  1  =yl  ,1  +yl,2  +  yl,2  +  y2,1  on  *  • 

The  differential  equations  (1.29)  may  suffer  a  loss  of  ellipticity 
at  solutions  that  have  sufficiently  severe  local  deformations.  The  el¬ 
lipticity  conditions  for  the  system  (1.29)  are  discussed  in  [7]  and  due 
to  their  complicated  nature  will  not  be  reproduced  here. 

The  special  compressible  material  that  concerns  us  has  the  elastic 
potential 


W(I,J)  =  |(IJ"2  +  2J  -4)  (v  >0)  ,  (1.33) 

where  u  is  a  constant.  This  strain  energy  density  was  proposed  by 
Blatz  and  Ko  [8]  to  model  a  highly  compressible  rubber-like  material. 

The  basic  properties  of  the  Blatz-Ko  material  are  investigated  in  [9]. 

It  is  shown  in  [9]  that,  for  the  material  characterized  by  (1.33),  the 
coordinate  equations  of  equilibrium  (1.29),  are  elliptic  at  a  solution 
ya  and  at  a  point  (x^Xg)  if  and  only  if  the  invariants  I  and  J, 
found  through  (1.32),  satisfy 

6  .  is  the  two-dimensional  Kronecker  delta. 
aB 


I 
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2  J  s  I  <  4  J 


(1.34) 


In  view  of  the  potential  loss  of  ellipticity,  it  is  natural  to  require 
of  the  coordinates  ya(x-|  only  relaxed  smoothness  specified  for 
u  on  &.  Thus,  there  may  be  curves  C  across  which  discontinuities  in 
the  deformation  gradients  y  occur.  Across  such  equilibrium  shocks, 

a  ,p 

the  Piola  tractions 


must  be  continuous  to  maintain  equilibrium. 

Finally,  the  homogeneous  deformation  corresponding  to  a  state  of 
plane  strain  uniaxial  stress  parallel  to  the  x2-axis  is  of  interest. 

For  the  Blatz-Ko  material  (1.33),  such  a  deformation  is  described  by 

y^r1^  ,y2=*x2  »  A>0  for  all  (x1,x2)€£>  ,  (1.36)1 


where  the  stretch  A  is  a  constant.  The  associated  true  and  nominal 
stresses  are 


t22  =  p(1-A  l  )  ,  Tn  =  Ti2  =  t21  =  ®  ^  ’ 

(1.37) 

-1/3  -3 

<?22=|i(A  —  A  )  ,  CT1 1  =  a12  =  a21  ~  ^  ^  • 

Suppose  a  solution  to  the  governing  differential  equations  (1.29) 
corresponds  locally  to  a  state  of  plane  strain  uniaxial  stress  with 
principal  stretch  ratio  A.  If  x  lies  on  the  interval  (Ag\A0), 
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where 

Xq=(2+^)3/4  , 

the  field  equations  will  be  locally  elliptic  (see  Fig 


t 
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2.  ANTI-PLANE  SHEAR  DEFORMATIONS  NEAR  THE  TIP  OF  A  CRACK 

In  the  present  section  we  give  the  formulation  and  asymptotic  solu¬ 
tion  of  a  nonlinear  crack  problem  involving  finite  anti -plane  shear  de¬ 
formations.  This  analysis  is  included  to  provide  valuable  insight  into 
the  more  challenging  tension  crack  problem  considered  in  the  next  section 
Let  the  open  cross-section  6  be  the  plane  domain  exterior  to  the 
line-segment 


£=  £ (xi  ^2) |-a  sx-|  ^a  ,  X£  =  0 }  ,  (2.1) 

as  depicted  in  Fig. 3.  Thus,  the  region  ft  corresponds  to  the  undeformed 
configuration  of  an  all-around  infinite  body  with  a  plane,  infinitely 
long  crack  of  width  2a. 

The  body  is  composed  of  the  special  incompressible  material  intro¬ 
duced  in  Section  1  and  characterized  by  the  response  curve  in  simple 
shear  (1.23).  Suppose  the  body  is  subjected  at  infinity  to  a  simple 
shear  parallel  to  the  edges  of  the  crack.  An  assumption,  consistent  with 
the  loading  and  the  particular  material,  is  that  the  deformation  is  en¬ 
tirely  one  of  anti-plane  shear  as  given  in  (1.14).  Referring  to  Eqs. 
(1.17),  (1.18)  and  (1.24)  we  obtain  the  governing  differential  equation 
for  the  unknown  out-of-plane  displacement  u: 

v^u  =  0  for  | vu |  <  1  ,  (2.2a) 

[|vu|"3/2u  q]  a  =  0  for  | vu |  >  1 


(2.2b) 
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At  infinity  we  stipulate 

u  =  yx2  +  o(1)  on  £  xaxa-oo  ,  (2.3) 

where  the  positive  constant  y  is  the  amount  of  shear.  Further,  we  as¬ 
sume  that  the  crack  boundary  in  the  deformed  body  remains  traction  free. 
According  to  (1.8),  (1.25)  and  the  inequality  in  (1.15),  this  require¬ 
ment  is  equivalent  to 


u>2(x1,0±)  =  0  (-a < x-|  <a)  .  (2.4) 

Note  that  the  solution  u  to  the  boundary  value  problem  (2.2)-(2.4) 
is  to  obey  the  continuity  requirements  set  down  in  Section  1.  Also,  for 
the  problem  at  hand  we  require  that  u  be  bounded  within  any  circle  of 
finite  radius  centered  at  the  crack-tips  and  the  limits  vu(x-|,0±)  are 
to  exist  and  be  continuous  for  -a<x-j  <a. 

The  nonlinear  crack  problem  formulated  above  is  one  of  considerable 
difficulty.  Knowles  and  Sternberg  [4]  considered  the  so-called  "small- 
scale  nonlinear  crack  problem"  associated  with  the  global  problem  de¬ 
scribed  above.  In  the  small  scale  problem,  y  is  assumed  to  be  small 
compared  to  unity  and  the  finite  crack  is  replaced  by  a  sem{- infinite 
one.  One  seeks  a  solution  of  (2.2)  which  satisfies  (2.4)  on  the  semi¬ 
infinite  crack  and  which,  at  infinity,  "matches"  the  near-tip  field  pre¬ 
dicted  by  the  linearized  theory  for  the  original  global  problem  associat¬ 
ed  with  the  finite  crack.  An  exact  solution  to  this  small-scale  prob¬ 
lem  was  constructed  in  [4]  by  means  of  a  hodograph  transformation.  Here, 
we  consider  the  global  problem  and  use  a  direct  asymptotic  approach  to 
study  the  field  near  the  crack-tip. 
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Let  (r,e)  be  local  polar  coordinates  in  the  undeformed  configura¬ 
tion  as  shown  in  Fig. 3.  Then, 

x-j-a  =  rcos0  ,  x2=rsine  .  (2.5) 

The  tip  of  the  crack  is  denoted  by  r=0,  while  0  =  it  and  0  = -rr  rep¬ 
resent  the  upper  and  lower  crack  faces,  respectively. 

We  first  investigate  the  possibility  of  "smooth"  solutions  near  the 
crack-tip.  Suppose  the  global  solution  to  the  crack  problem  admits  the 
asymptotic  representation 


u~rmv(0)  as  r-0  ,  (2.6)^ 

uniformly  on  -irsesi.  In  (2.6)  m  is  a  constant  in  the  range 

(Um<l  ,  (2.7) 

and  v  is  a  "smooth"  function  on  [-ir,ir].  Here,  we  say  a  function  is 
smooth  if  it  is  at  least  twice  continuously  differentiable  on  its  domain 
of  definition.  In  the  present  problem,  symmetry  implies 

v(e)  =  -v(-0)  (-ir£0sSTr)  .  (2.8) 

The  boundary  condition  (2.4)  yields 

v  =  0  on  0  =  -it, it  ,  (2.9) 

where  the  dot  denotes  the  differentiation  with  respect  to  0.  The  re¬ 
strictions  (2.7)  on  m  guarantee  that  the  displacement  remains  bounded 

Hhe  asymptotic  equality  symbol  is  used  in  the  following  sense; 
u~rmv(e)  is  equivalent  to  u'^vtol+ofr1’1)  as  r-0. 
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near  r  =  0  for  -irsesir,  while  permitting  unbounded  displacement  gra¬ 
dients. 

Equations  (2.6)  and  (1.18)  imply 

|vu|  ~p(0)rm"1  as  r-0  ,  (2.10) 

where 

P(e)=/m2v2(e)  +  v2(e)1  ,  (2.11) 

for  -irseiir.  Clearly,  if  p(e)>0,  then  |vu|  >1  as  r-0,  -irsesir, 
so  that  (2.2b)  applies.  Assume  this  to  be  the  case  and  substitute  (2.6) 
into  (2.2b).  Retaining  only  the  dominant  power  of  r  provides  the  gov¬ 
erning  differential  equation  for  v(e): 

2pv  -  3pv +  m(3  -  m)pv  *  0  ,  (2.12) 


on  (-ir, it)  . 

If  m=0  in  (2.12),  we  obtain  the  solution  v(e)  =  ce  +  b  (b  and  c 
are  constants).  Consideration  of  (2.8)  and  (2.9)  requires  v(e)*0  for 
-irsesir.  We  may  thus  restrict  m  to  the  interval  (0,1). 

The  differential  equation  (2.12)  can  be  analyzed  in  the  phase  plane. 
Details  of  this  analysis  are  given  in  [12].  One  finds  in  this  way  that 
the  general  solution  of  (2.12)  may  be  written  in  the  form 

mv(e)  =  V cos  *(e)U0  +  cos  2if»(e) |'1/2(1  _m)  »  (2.13) 

where  ij>(e)  is  the  general  solution  of 

(1/3  +  COS  2tp) ip  +  (aiQ  +  COS  2i|»)  *  0  , 


(2.14) 
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for  -tt  < ©  < 7r .  In  (2.13),  V  is  constant  which  we  assume  to  be  positive 
and 


w  »  1/3(3  -2m)  .  (2.15) 

Observe  that  for  0<m<l,  one  has  1/3<uq<1.  In  Eq . (2.13)  we  must 
have  uo  +  cos  2<p  f  0  to  assure  bounded  displacements.  Further,  we  note 
that  if  p  is  calculated  from  (2.13),  (2.11),  it  vanishes  nowhere  on 
[-ir.ir],  as  was  assumed. 

In  view  of  Eqs.  (2.8),  (2.13)  and  (2.14)  we  choose 

tf(0)  =  */2  ,  (2.16) 

so  that  (2.9)  holds  if  and  only  if 

4>(ir)  =  2nir  ,  ij)(-ir)=(2n  +  l)ir  ,  n  =  0,±1  ,±2, . . .  (2.17) 

In  Fig. 4  we  sketch  the  curve  ip  =  4#(e )  governed  by  (2.14),  (2.16)  and 
(2.17)  for  the  case  n=0.  It  is  apparent  from  the  figure  that  there 
are  no  continuous  solutions  of  the  boundary  value  problem  for  \p  -  this 
result  is  found  to  be  independent  of  the  choice  of  n.  Consequently, 
there  does  not  exist  a  smooth  function  v  on  [-ir.Tr]  that  satisfies 
(2.12),  (2.8)  and  (2.9). 

We  now  wish  to  generalize  the  Ansatz  (2.6)  concerning  the  form  of 
u  near  r  =  0  in  such  a  way  as  to  permit  discontinuities  in  vu  across 
certain  curves  issuing  from  the  crack-tips.  Accordingly,  we  suppose 
that  there  are  two  curves,  S+  and  S",  defined  in  a  neighborhood  of 
the  origin  by 


L. 


I 
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8 ±:  e*±e(r)  ,  0<rirQ  ,  (2.18) 

where  e(r)  is  a  smooth  positive  function  such  that 

e(r)~  Ars  as  r-0  ,  (2.19) 

and  A>0,  s^O  are  constants  to  be  determined.  If  s  =  0  in  (2.19), 
then  A  is  required  to  satisfy  0<A<ir  and  the  curves  S*  are  tangent 
to  the  rays  e  =  +A  at  the  origin.  When  s>0,  S+  and  S"  are  both 
tangent  to  the  x^axis  at  the  origin.  Since  our  interest  lies  in  re¬ 
producing  the  results  in  [4],  in  which  the  shocks  Sr  are  found  to  be 
tangent  to  the  x^-axis  at  r=0,  we  treat  here  only  the  case  s> 0  in 
(2.19).  The  case  s  =  0  is  discussed  in  Appendix  A. 

It  is  convenient  to  introduce  the  regions  AC,  £  and  C  (Fig. 5) 
as  follows 

Ac  =  {(r,e)|-e(r)  <e  <e(r)  ,  0<r£rQ3  , 

5+=  ((r,9)ie(r)  <9<iv  ,  0<rsrQ]  ,  >.  (2.20) 

c"- {{r,e)|-ir<e<-0(r)  ,  0<rsrQ3  . 

We  first  investigate  solutions  of  the  differential  equation  (2.2b)  which, 

in  the  region  AC,  have  the  asymptotic  form  (2.6),  where  v(e)  is  a 

smooth  odd  function  defined  for  -e.sese,  for  some  e  .  0<e  <n. 

oo  oo 

The  exponent  m  is  now  permitted  to  be  neaative,  but  rmv(e)  must  re¬ 
main  bounded  as  (r,e)  approaches  the  origin  from  within  fc. 

If  m<0,  we  find  that  (2.10),  (2.11)  still  apply  as  r-0  in  * 
and  v(e)  must  again  satisfy  (2.12),  which  leads  to  the  same  implicit 
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r 


representation  given  in  (2.13)-(2.15)  except  that  now  m<0.  Note  that, 
in  this  case,  wQ>l  and  hence  <uQ  +  cos  2*  f  0.  As  before,  (2.16)  ac¬ 
counts  for  the  parity  condition  (2.8).  Integrating  (2.14),  subject  to 
(2.16),  yields  <p(0)  implicitly  through 


6  a  -(♦  -  */2)  +  jjrJ  w-jtan"1  [c^tan^  -  */2)] 


(“0o<e<eo)  ,  (2.21 ) ^ 


where 

.ra-T7*  , 

r[-T ]  <m<°> 


and 


(2.22) 


e0  =  -1/2 cos'1  (|)  +  tan'1  ,  0<eQ<ff  .  (2.23) 

At  this  point  Eq . (2 . 21 )  cannot  be  inverted  explicitly  to  furnish 
*a*(0)-  rt  does.  however,  indicate  that  *(e)  is  a  continuously  dif¬ 
ferentiable  function  which  increases  monotonically  on  (-0  ,0  )  from  a 
value  </»0=l/2cos  ^  (-1/3)  to  it  -  ipQ  and  is  antisymmetric  about  e  =  0, 
*  =  ir/2.  In  view  of  these  facts,  we  confirm  through  (2.13),  (2.21)  the 
existence  of  a  smooth  function  v  on  (-0Q,eo)  which  satisfies  (2.12) 
and  has  the  appropriate  symmetry.  Thus, 

u~rmv(e)  (m<0)  as  r-*0  ,  (2.24) 

represents  an  asymptotic  solution  to  the  hyperbol ic  differential  equa¬ 
tion  (2.2b)  on  V. 

’Here,  the  inverse  trigonometric  functions  take  their  principal  values. 


a‘— '■ 
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For  |e(«1,  Eqs.  (2.13)  and  (2.21)  provide 

v  =  pQ C 8  +  g- ( 3  -  4m ) 6 3  0(65)]  as  6-0  ,  (2.25) 

where 

Po  =  p(0)  .  (2.26) 

From  (2.18),  (2.19),  (2.24)  and  (2.25)  we  observe  that  the  displacement 
remains  bounded  as  the  crack-tip  is  approached  from  within  V  if 
s  +  m^O. 

We  now  seek  a  solution  to  (2.2)  valid  on  the  region  £+.  Assume 
the  displacement  on  £+  admits  the  general  form  (2.6)  where  the  expo¬ 
nent  miO  takes  on  the  smallest  possible  value  -  boundedness  of  the 
displacement  as  the  origin  is  approached  in  £+  necessitates  m  be 
nonnegative.  We  find  that 

u~b  +  rv(e)  as  r-0  on  s+  ,  (2.27) 

where  b  is  a  positive  constant  and  now  v  is  a  smooth  function  on 
(0,w]  that  obeys  the  boundary  condition  v(n)=  0.  Note  that  the  con¬ 
stant  b  is  the  solution  corresponding  to  m=0  and  reca  ,  for 
0<m<l  in  (2.6),  there  are  no  smooth  solutions  of  (2.2)  on  (0,tt]  sat¬ 
isfying  the  boundary  condition  on  e  =  ir. 

Substituting  (2.27)  into  either  (2.2a)  or  (2.2b)  leads  to  the  same 
differential  equation: 


v  +  v  =  0  (0<  6<  it)  . 


(2.28) 


The  most  general  solution  of  (2.28)  for  which  the  boundary  condition 
v(e) =  0  is  satisfied  is 

v(e)  =  ccose  ,  0<esit  ,  (2.29) 

where  c  is  a  constant.  In  C,  the  dominant  term  analogous  to  (2.29) 
is  obtained  by  symmetry. 

It  is  convenient  at  this  point  to  present  the  results  of  a  higher 
order  analysis.  If  | c |  < 1  we  find 

u  ~b  +  cr  cos  e  +  dr  cos  k(e  -  rr)  (k>l)  as  r-0  on  5  .  (2.30) 

Note  that  | vu |  — | c |  < 1  as  r-*0  which  assures  the  ellipticity  of  the 
displacement  equations  of  equilibirum  on  £+. 

What  remains  now  is  to  construct  an  asymptotic  solution  to  the  crack 
problem  near  the  crack-tip  -  continuous,  piecewise  smooth  and  with  con¬ 
tinuous  nominal  tractions  across  S±.  To  this  end,  we  confine  our  atten¬ 
tion  to  0^0  i it  and  note  that  (2.8)  generates  the  solution  in  the  lower 
half-plane. 

The  requirement  that  u  be  continuous  across  S+,  together  with 
(2.18),  (2.19),  (2.24),  (2.25)  and  (2.30),  clearly  implies 

b~ApQrfrHn  ,  (2.31) 

so  that  necessarily 

A=p0^b  ,  8  =  -m  (2.32) 

We  now  address  the  continuity  of  the  non-zero  nominal  traction  s^ 


-22- 


(1.25)  across  the  shock  S+.  The  asymptotic  forms  of  the  traction  (1.25) 
on  the  two  sides  of  S+  are  evaluated  through  (1.18),  (1.24)  where  the 
displacement  on  the  hyperbolic  side  of  the  shock  is  given  by  (2.24), 

(2.25)  while,  on  the  elliptic  side,  we  use  (2.30).  The  tractions  are 
found  to  be  continuous  across  S+  if  and  only  if 

Py2rl/2(m-l)^(m_  i)bcr‘l +dPoksin  k7rrk+ni"2  .  (2.33) 

The  possible  relationships  between  the  exponents  k  and  m  that 
permit  an  asymptotic  balance  in  Eq.(2.33)  are  described  in  Fig. 6.  Asymp¬ 
totic  solutions  to  the  next  order  of  the  displacement  equation  of  equi¬ 
librium  (2.2a),  (2.2b)  were  found  on  £+  and  X  in  an  attempt  toelim- 
inate  some  of  the  cases  I-V  in  the  figure.  However,  the  solutions  had 
sufficient  flexibility  so  as  to  allow  the  necessary  higher  order  match¬ 
ing  across  an  appropriately  chosen  shock  S+  for  all  the  five  cases. 

Our  purpose,  in  this  paper,  is  to  reproduce  the  solution  in  [4]  in  the 
vicinity  of  the  crack-tip  and  hence  we  assume  case  II  to  hold,  so  that 

m=  -1  ,  k  =  2  .  (2.34) 

Equations  (2.33),  (2.34)  reveal  that 

/  PV2\ 

-l<[c=--|r-]<l  .  (2.35) 

We  note  that  for  m=-l,  (2.21)  can  be  inverted  analytically  and  on  sub¬ 
stitution  into  (2.13)  furnishes,  after  some  algebra 
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v  ( e )  =  2 


,2-5/2  (3cose*P(9))3/f 


(cos  e  +  P(e) ) 


T7? 


■sin  9 


-9«<9<9J 

O  0 


(2.36) 


where 


P(e)  =  (9cos2  e- 8)1/2  ('e0<9<90)  »  (3.37) 

and  eQ =  cos"^ (2/2/3).  This  completes  the  matching  to  leading  order. 

To  conclude,  we  summarize  the  asymptotic  representation  for  the 
solution  to  the  crack  problem  (2.2)-(2.4)  in  the  vicinity  of  the  right 
crack-tip.  Equations  (2.24),  (2.30)  together  with  (2.34),  (2.35)  give 


1/2 

PQ  2  +  ^ 

u~b  -  r  cos  6  +  dr  cos  2e  as  r-0  on  £  , 


2b 


u~r_1v(e)  as  r-0  on  V  , 


(2.38) 


where  b>0,  0<po<«/2F  and  d  are  unknown  constants  and  v(e)  is  given  in 
(2.36).  In  (2.38),  C  is  the  elliptic  region  and  K  the  hyperbolic 
region  described  in  (2.20).  The  displacements  on  s"  are  obtained 
through  symmetry.  Equations  (2.18),  (2.19)  and  (2.32)  indicate  the  elas- 
tostatic  shocks  S+  and  S"  have  the  following  asymptotic  form 


S": 


as 


r-0  . 


(2.39) 


The  above  results  are  consistent  with  those  found  in  [4]  on  taking 
the  limit  r-0.  However,  in  the  local  analysis  presented  here,  the 
constants  b,  d  and  pQ  in  (2.38)  remain  undetermined.  In  contrast, 
the  specific  small-scale  nonlinear  crack  problem  treated  in  [4]  leads 
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to  fully  determined  values  of  b,  d  and  pQ.  Furthermore,  for  each  suc¬ 
cessively  higher  order  calculation  we  obtain  one  additional  unknown  con¬ 
stant.  The  reader  may  refer  to  [4]  to  obtain  representations  for  the 
stresses  near  the  crack-tip. 
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3.  PLANE  STRAIN  DEFORMATIONS  NEAR  THE  TIP  OF  A  CRACK 

In  this  section  we  present  the  formulation  and  first  order  asymp¬ 
totic  analysis  of  the  nonlinear  plane  strain  crack  problem  which  consti¬ 
tutes  the  main  purpose  of  this  study. 

Let  ft  be  the  undeformed  configuration  of  a  cylindrical  body  whose 
cross-section  &  in  the  (x-|  ,x2)-plane  is  described  in  Fig. 3.  The  line 
segment  X  (2.1)  is  the  boundary  of  X  and  represents  a  crack  of  width 
2a.  The  material  making  up  the  body  is  compressible  and  has  the  Blatz-Ko 
plane  strain  elastic  potential  (1.33). 

Suppose  the  body  is  subjected  to  uniform  uniaxial  tension  at  infin¬ 
ity  perpendicular  to  the  crack  faces,  so  that  the  actual  stresses  satis¬ 
fy 

TaB(xl’x2)  =  <s2a,V  +  o(1)  as  VcT00  *  (3J) 

where  t>0  is  the  magnitude  of  the  applied  load.  The  crack  faces  are 
to  remain  traction  free.  Applying  (1.7)  and  (1.8)  results  in  the  bound¬ 
ary  conditions 


CTa2 (X1  »®± )  =  0  (_a<xi<a)  •  (3.2) 

As  a  consequence  of  (3.1)  and  (3.2)  the  deformation  conforms  to  (1.26), 
corresponding  to  plane  strain  parallel  to  the  (x-|  ,x2)-plane.  Accordingly, 
the  relevant  coordinates  ya  must  satisfy  the  equations  of  equilibrium 
(1 .29),  together  with  (1.32)  and  (1.33).  The  conditions  at  infinity  (3.1 ) 
and  boundary  condition  (3.2)  are  represented  in  terms  of  the  coordinates 


through  (1 . 30 ) - ( 1 .32).  The  solution  y  to  the  boundary  value  problem 

a 

must  be  continuous  and  piecewise  smooth  on  6  as  indicated  in  Section  1. 
We  insist  also,  that  the  solution  be  bounded  in  the  vicinity  of  the 
crack-tip  and  the  limits  y  .  (x1t0±)  exist  and  are  continuous  for 

01 »  P  I 

-a<x1  <a. 

An  analytical  solution  to  the  global  problem  is  not  attempted  here. 
We  consider  instead,  the  asymptotic  character  of  the  solution  in  the 
neighborhood  of  the  crack-tipJ 

Let  (r,e)  be  the  polar  coordinates  introduced  in  (2.5)  and 
ya  =  yct(r,e)  be  the  local  deformation  near  the  right  crack-tip.  From 

O 

(1.29),  with  W  given  by  (1.33),  we  obtain  the  appropriate  form  for  the 
governing  differential  equations,  namely 

A-2<Va.r*r'2V<M) 

+  (-)ar‘1J3(Rjeye>r-R>ryB>6)  =  0  (a^e)  ,  (3.3) 

for  r>0,  — ir < e < ir.  In  (3.3), 

R(r,e)  =  1  -  IJ"3  ,  (3.4) 


where  (1.32)  gives 


O(r.0).r-1(y,iry29-yli9y2jr)>O  ,  (3.5) 

I(r,6) -yf>r  +  y|>r*r-2(y^e  +  y|je)  ,  (3.6) 


T 


The  existence  of  a  solution  is  assumed. 
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for  r>0,  -irs0sii.  The  loading  at  infinity  (3.1)  indicates  the  fol¬ 
lowing  coordinate  symmetries: 


y^r.e)  =  y1(r,-e)  ,  y2(r,e)  =  -y2(r,-e)  (r>0,  -Tr^e**)  .  (3.7) 

Substituting  (1.31),  (1.33)  into  (3.2)  provides  the  boundary  conditions 
on  the  crack  faces 


n.rJ'l.o’^ya.e’0 


Von  6  =  -it, it  for  r>0  .  (3.8) 


To  complete  the  formulation  we  stipulate  that  the  resulting  deformations 
have  the  same  smoothness  as  that  specified  for  the  global  problem. 

Now  turn  to  the  asymptotic  solution  of  the  crack  problem.  We  fol¬ 
low  the  anti-plane  shear  example  in  the  previous  section  and  assume  that 
the  local  deformation  field  is  represented  by 


m 

ya~r  ava(9)  as  r-*0  (no  sum  on  a)  ,  (3.9) 

where  ma  are  constants.  In  Aooendix  B  we  consider  0  smin[m^  ,m2] <  1 
in  an  attempt  to  find  a  solution,  consistent  with  (3.3)-(3.8),  in  which 
va  are  smooth  functions  on  [-*,ir].  We  conclude  from  the  analysis  that 
no  such  deformations  exist.  Motivated  by  the  anti-plane  shear  problem, 
we  now  let  8+  and  8“  be  the  curves  originating  at  r  =  0  described 
in  (2.18),  (2.19).  Also,  define  the  regions  U,  and  e"»  in  the 
neighborhood  of  the  crack-tip  through  (2.20)  (see  Fig. 5). 


Suppose  (3.9)  holds  on  U  with 
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-  m<  0  ,  m<m^  ,  (3.10)^ 

and  va  are  smooth  functions  on  (-90»90)  (0<eo<ir).  We  assume  here, 

that  the  deformations  are  "most  severe"  in  the  Xg-direction.  We  let 

J~rvq(e)  as  r-0  on  V  ,  (3.11) 

where  q  is  a  continuously  differentiable  function  for  -e0<e<e0  and 
(3.5),  (3.9)  and  (3.10)  infer  the  constant  v>2(m-l).  Further,  (3.6), 
(3.9)  and  (3.10)  provide  the  asymptotic  form  of  the  scalar  invariant  I, 
namely 


I~r2(m-1)p2(0)  as  r-0  on  V  ,  (3.12) 


where 

p(e)  =ym2v2(e)  +  v^e)  on  (-eo,eo)  .  (3.13) 

Proceeding  as  in  the  case  (B.23)  in  Appendix  B  we  find  that  (B.25)- 

(B.27)  hold  on  (-0O»9O)  for  m<0.  Equations  (B.29)-(B.32)  provide 

expressions  for  Vg  and  q  on  ("90»90)  while  the  symmetry  (3.7)  and 

smoothness  across  the  ray  9  =  0  is  guaranteed  by  (B.34).  In  (B.34)  we 

2 

choose  e-|  =  0  and  hence  (B.29)-(B.32)  reduce  to 

mvg(e)  =  Vgcos  i)(e)|c0 +  cos  2^(e)|^2^m“^  ,  (3.14) 

Hhe  case  m1sm2<0  is  treated  in  a  manner  similar  to  that  of  the  case 
of  positive  equal  exponents  discussed  in  Appendix  B.  We  find  that  the 
coordinates  (3.9)  and  the  parity  condition  (3.7)  are  incompatible. 

2 

The  case  ei^l  Is  not  treated  here.  We  find  this  leads  to  inconsis¬ 
tent  result  J(r,0)  =  0  for  aJJ_  sufficiently  small  r>0. 
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q(e)  =  Q|e0  +  COS2*(e)|1/3(m_1)  ,  (3.15) 

for  -0O<0<0O>  where  V2  and  Q  are  positive  constants  and  ^ ( 9 ) 
satisfies 

(1/2  +  cos  +  {eq  +  cos  2\(i)  =  0  ,  ip(0)  =  tt/2  .  (3,16) 

Furthermore, 

eq  =  1  -  m/2  ,  v-  2/3(m-  1)  .  (3.17) 

Integrating  (3.16)  yields 

0  = -(</;  -  tt/2) +||-5-^j-E4tan_1  [e4tan(i<)  -  ir/2)]  (-0o<0<eo)  .(3.18) 

where 

1  /2 

C4'!1—]  !rn<0)  ,  (3.19) 

and 

9o  =  -ir/6  +  {^m}E4tan  ^  ( e4tan  tt/6 )  ,  0<eo<ir  .  (3.20) 

The  expressions  (3.18)-(3.20)  indicate  that  1(1(0)  is  a  smooth  function 
on  ("e0*e0)  and  hence  v2(e)  (3.14),  q(o)  (3.15)  are  also  smooth 

with  the  correct  parity  for  -0O<0<0O  • 

We  now  evaluate  the  y^ -coordinate  to  leading  order.  Note  that 
(3.13),  (3.14)  and  (3.15)  imply 

q(0)  =  QoP2/3(0)>O  on  (-0O,0Q)  , 


(3.21) 


with  Qo  =  QV22/3.  The  Jacobian  is  evaluated  using  (3.5),  (3.9),  (3.10) 
and  again  through  (3.11),  (3.21).  Equating  the  two  results  provides  the 
asymptotic  equality 

r^1  (m1v1v2-mv1v2)~Q0r2/3^m_1^p2/3  as  r-0  ,  (3.22) 

m,  /m 

on  V.  Assume  first,  m1<l/3(4-m)  and  find  v1  (e )  =  V,  |  v2(e)| 
on  ( -6q , 6q )  where  V^(>0)  is  a  constant.  The  restrictions  (3.10) 
imply  that  m-j  =  0  and  v-|(e)  =  Vi  is  the  only  smooth  solution  across 
e  =  0.  Taking  y-j(0,e)  =  0  (-ir£6*ir)  requires  V-j  =  0 ?  Now  let 

m1  =  1/3(4  -m)  ,  (3.23) 


and  obtain  a  non homogeneous  first  order  differential  equation  for  v1 
from  (3.22).  We  establish  that  the  homogeneous  solution  is  unbounded  on 
e  =  0  and  is  neglected,  while  the  smooth  particular  solution  gives 


»(e)  | 


let 


-Ec 


p2/3(?)vo  3(<p)d9  on 


(-eo,eQ)  ,  (3.24) 


where 

e5 s  <  1  for  m<0  •  (3.25) 

Equations  (3.9),  (3.10),  (3.14),  (3.17),  (3.23)-(3.25)  provide  an 
asymptotic  solution  to  the  differential  equation  (3.3)-(3.6)  such  that 


T 


The  deformation  is  normalized  so  that  the  crack-tip  does  not  move  with 
respect  to  the  xi-axis. 
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Mm<0)  as  r-0  on  V 


(3.26) 


y2  ■  h2~rmv2(e)  J 

Further,  (3.11),  (3.17)  and  (3.21)  give 

OV^-V/^)  as  r-0  on  V  ,  (3.27) 

while  substituting  (3.26)  and  (3.27)  into  (3.3)  (a  =  l)  reveals 

R=  1  -  Q‘3  +  0(r8/3^1-rn^)  as  r-0  on  U  .  (3.28) 

It  is  apparent  from  (3.12)  and  (3.27)  that  I»J  on  V  and  thus  the 
coordinate  equations  of  equilibrium  are  hyperbolic  at  the  solution  (3.26) 
on  U. 

For  |e|«l,  (3. 14)-(3. 16)  and  (3.24)  provide 


3q„p;1/3.  , 


.2  .  _/ -T'' 


vr  (4-°m)  tl-g»(4-»)e‘»0(»’)] 


v2  =  pQ[e  -  ^-m(5m  -  4 ) e3  +  0(e5)] 


v  as  0-0  ,  (3.29) 


q  =QqPq/3[1 -ym(m-l)e2+O(04)] 


J 


where 

PQ  =  P(0)  .  (3.30) 

Equations  (2.18) ,  (2.19),  (3.26)  and  (3.29)  infer  that  the  y2-coordinate 
remains  bounded  on  U  in  the  limit  r-0  if  g  +  maO. 
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We  now  seek  solutions  to  the  governing  equations  (3.3)  on  £+.  Sup- 


y  =t  ~b  +  rv  (e)  as  r-  0  on  s+  , 
a  a  a  a 


(3.31) 


where  b  are  constants  and  v  are  now  smooth  functions  on  (O.irl.  As 
a  a 

in  the  anti-plane  shear  problem  we  obtain  the  representation  (3.31)  after 
first  assuming  (3.9)  holds  on  £+  with  0 smin[m^  .ir^] <  1  and  draw  on 
results  from  Appendix  B.  The  constants  b&  are  the  solutions  (B.5) 
for  m^=m2  =  0.  The  case  m.|  =  =  m  (0<m<1)  yields  v-|(0)  =  cv2(e) 

on  (0,tt] ,  where  c  is  a  constant.  If  c  =  0  we  recover  the  unequal 
exponent  case.  The  boundary  conditions  on  e  =  it  are  given  in  (B.28). 
After  a  tedious  analysis  we  prove  that  no  smooth  solution  v2  exists 
on  (0,ir]  that  is  compatible  with  (B.28)  and  can  be  matched  to  the  hyper¬ 
bolic  solution  on  V  across  8+. 

Substituting  (3.31)  into  (3.3)-(3.6)  and  retaining  the  leading  order 
terms  gives 


[(v^  +  2v2)v2  + v]v2v1](v1  +  V'j ) 

-C(v2  +  2v^)vi  +v1v2v2](v2  +  v£)  =  0  , 

[3(v*  +  v^Vg  +  (v-j V-|  +  v2v2)v2](v1  +  v] ) 

-[3(v^  + Vg)v1  +  (v]v1  +  v2v2)v1](v2  +  v2)  =  0 


(3.32) 


for  0 < 0 < it.  The  differential  equations  (3.32)  are  equivalent  to 


v  +v  =0  on  (0,ir)  , 

a  a 


(3.33) 


i 


unless  the  determinant  of  the  coefficients  of  (v  +v  ) 

a  a 

(0,ir)J  Assume  (3.33)  holds  and  thus 

Vj  (e)  =  x.|Cos  e  -  x2(x3  +x  2)"2/,3sin  e  , 

2  2  -2/3 

v2(e)  =  x2cos  e  +  x1(x1 +x  2)  sine  , 


vanishes  on 


>  (3.34) 


for  0<6iir,  where  x  are  constants.  The  solution  (3.34)  obeys  the 

at 

boundary  conditions  (3.8). 

Observe  from  (3.26),  (3.29),  (3.31),  (3.34)  and  (2.18)  that 
h-j  ~3^o*>o^^4  “  m)-1^3^  "  ,  a^~b-|  +  X^r  as  r-0  on  8+.  As 

l/3(4-m)>l  for  m<0  we  arrive  at 

b|  =  0  ,  Xj  =  0  ,  (3.35) 

as  a  prerequisite  for  continuous  coordinates  across  S+.  Letting  the 
scalar  invariants  take  on  their  leading  order  values  calculated  through 
(3.5),  (3.6),  (3.31),  (3.34)  and  (3.35),  the  ellipticity  condition  (1.32) 
reduces  to  (2  -  ^J)3^<  |x2|  <  (2  + /I)3^.  Physically,  in  the  vicinity 
of  the  crack  face,  we  expect  the  leading  order  homogeneous  deformation 
(3.31),  (3.34)  and  (3.35)  to  represent  a  state  of  uniaxial  tension  par¬ 
allel  to  the  y2~direction.  This  limits  the  range  of  the  stretch  ratio 
further,  so  that  elliptic  solutions  of  interest  obey 

1  <X<  (2+ /J)3/4  ,  x  =  -x2>0  .  (3.36) 


T 


It  can  be  shown  that  functions  v0  that  produce  a  vanishing  determi¬ 
nant  are  not  consistent  with  either  the  boundary  conditions  on  the 
crack  face  or  the  matching  across  S+. 


-34- 


In  what  follows  we  require  the  solution  on  £+  to  higher  order. 
Noting  (3.31)  and  (3.34)-(3.36),  we  now  assume 


=  e.j  ~x"^3r  sin  e  +  r^f^  (e) 
4*  k 

y2  =  e2  -  Xr  C0S  0  +  r  '2^ 


as  r-*0  on  C  ,  (3.37) 


where  k>l  is  constant  and  f  are  smooth  functions  on  (0,ir]. 

a 

Equation  (3.37)  satisfies  the  differential  equations  (3.3)  to  the  appro¬ 
priate  order  if 


La(fl  ’f2;k,X)  =  0  on  » 


(3.38) 


where  La  is  the  differential  operator, 

L]  (^  ,f2;k,x)  ■  (2  +  K+  Kcos  29)^  +  2K(k  -  1  )sin  29^ 

+  k[k(2  +  K)  -  K(k  -  2)cos  2e]f -j 

+  N[sin  2ef2  -  2(k  -  1  )cos  2ef2  -  k(k  -  2)sin  2ef2]  ,  (3.39) 


l2 ( fi  ,f 2 ; k , A )  =  (2  +  M  -  M cos  2e)f2  -  2M(k  -  1  )sin  2ef2 
+  k[k(2  +  M)  +  M(k  -  2)cos  2e]fz 

+  N[sin  29^  -  2(k  -  1  )cos  29^  -  k(k  -  2)sin  2ef^]  .  (3.40) 


In  (3.39)  and  (3.40) 


K  *  3X8/3-  1 


M=  3X‘8/3-  1 


=  X4/3  +  X_4/3 


(3.41) 
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The  boundary  conditions  (3.8)  infer 

3f1(1r)-X_4/3kf2(ir)  =  0  ,  f2(1r)-X_4/3kf1(1r)  =  0  .  (3.42) 

Furthermore,  (3.7)  implies 

y-j  =  i-j  =  ^  (r,-  e)  ,  y2  =  e2  = -e2(r,- e)  on  C  .  (3.43) 

We  now  generate  the  first  order  asymptotic  solution  with  the  appro¬ 
priate  smoothness  to  the  small  scale  nonlinear  crack  problem  (3.3)-(3.8). 
Consider  O^e^iT  and  refer  to  (3.43)  to  provide  the  deformation  in  the 
lower  half-plane.  Accordingly,  the  matching  conditions  are 

e  =  h  on  8+  ,  (3.44) 

a  a 

together  with  continuous  Piola  tractions  across  S+.  We  compute  the 
general  expression  for  the  nominal  tractions  through  (1.7),  (1.31)  and 
(1.33): 


V“'J'2(>,a.rnr+r'1 


ya.en8)*(-)C,R(yB.rne-r',)'6,9nr) 


(3.45) 


(n  .n„)  are  the  components  of  the  normal  vector  n  in  polar  coordinates. 

i  v 

Let  a  =  2  in  (3.44)  and  use  (2.18),  (2.19),  (3.26),  (3.29)  and 
(3.37)  to  find 


b2<~Apors+m  as  r-  0  . 


(3.46) 


Take  b2>0,  then 


^2  —  Ap0  9  8  ~  -m 


(3.47) 
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Similarly,  let  ct  =  1  in  (3.44)  and  use  (3.47)  to  obtain 

r]/V"  +  f]  (0)rk-1  ~  r'1/3(m_1)  as  r-0  .  (3.48) 

Here  and  in  what  follows  we  assume  that  the  functions  f  are  analytic 

a 

on  [0,ir]  and  hence  can  be  represented  by  a  Taylor  series  about  e  =  0. 
Before  matching  tractions  on  S+  we  modify  the  result  (3.28)  so 

that 


Q0  =  1  and  R  =  0(r8/3^1_mb  as  r-0  on  k  .  (3.49) 


We  confirm  (3.49)  by  first  assuming  the  existence  of  an  equilibrium  so¬ 
lution  in  the  neighborhood  of  the  crack-tip  that  has  the  general  form 
(3.26)  and  (3.37)  with  the  curve  S+  given  by  (2.18)  and  (3.47).  Force 
equilibrium  in  the  x-j-direction  on  a  circular  region  in  £>  centered  at 
r  =  0  necessitates  that  (3.49)  holds  on  taking  the  limit  r-oJ 

Equating  the  tractions  (3.45)  across  S+  and  drawing  on  (3.47)  and 
(3.49)  yields 


3^(0)-  x'4/3kf2(0)  =  0  , 


A’1/3(m-  1)(a‘8/3-  l)Ar'm 

+  x_4/3(f2(0)  -  A_4/3kf1(0))rk‘1  ^p’1/3r"1/3'm"1) 

The  relationships  between  the  exponents  k  and  m 


>  (3.50) 


as  r-0  .J 
that  could  provide 


— 

The  result  (3.49)  assures  no  concentrated  forces  act  at  the  crack-tip. 
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a  solution  to  (3.48)  and  (3.50)  are  presented  in  Fig. 7.  We  note  the 
similarity  between  this  figure  and  the  diagram  obtained  for  the  anti- 
plane  shear  problem  (see  Fig. 6).  Again,  we  choose  k  and  m  through 
a  distinguished  limit  and  thus 

m=  -1/2  ,  k  =  3/2  .  (3.51) 

In  view  of  (3.51),  the  equations  (3.48)  and  (3.50)  give 

A  =  xV3(f  po1/3  "  fl (0) )  *  (3.52) 

f--,  (o)  -  ^x'4/3f2(o)  =  o  ,  f2(o)-fA4/3f1(o)  =  r4/3P;;1/3 .  (3.53) 

Recall  in  (2.18),  (2.19)  that  A  is  a  positive  constant  and  hence 

fi(0)<!po1/3  •  (3-54) 

Equations  ( 3 . 38 ) - { 3 . 42 ) ,  (3.51)  and  (3.53)  constitute  a  boundary 
value  problem  for  f  on  [0,*].  For  convenience  set 

T  (0)  =  Py3fa(e)  (OSes tt)  ,  (3.55) 

where  f  satisfies 
a 

La(fl,f2;3/2,A)  =  0  on  (O.tt)  ,  (3.56) 


with 
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f1(o)-yx“4/3f2(o)-o  ,  f2(o)-|x4/3f1(o)  =  x'4/3 , 

(u)  -  jX"4//3f2(ir)  =  0  ,  f^Tt)  -  fx-4^3^  (it)  =  0  . 


(3.57) 


A  closed  form  solution  to  the  differential  equation  (3.56)  appears  to  be 
difficult  to  construct.  However,  a  numerical  analysis  indicates  that 
for  values  of  X  obeying  the  inequality  (3.36)  the  nondimensional  bound¬ 
ary  value  problem  (3.56),  (3.57)  does  possess  analytic  solutions.  For 
x  f  x*  ,  where 


X*=  2.519 


(3.58) 


the  solutions  depend  uniquely  on  X.  We  find  that  X =  x*  is  an  eigen¬ 
value  associated  with  the  homogeneous  version  of  (3.56),  (3.57).  Only 
one  eigenfunction  exists  together  with  the  particular  solution.  We  note 
that  in  all  the  solutions,  f^(0)<2/3,  complying  with  (3.54) . 1 

We  now  present  a  surmary  of  the  asymptotic  solution  to  the  crack 
problem  considered  here.  Applying  the  results  (3.47),  (3.51)  and  (3.53) 
to  (3.26),  (3.37)  and  (3.43)  gives 


^  ~x"1/3r  sin  e  +  p^V^  (e) 

-  *r  cos  e  +  p"^r3/^f„(e) 


>  as  r-0  on  £*  ,  (3.59) 


^Analytic  solutions  to  (3. 38)-(3 .42)  and  (3.53)  were  found  for  integer 
values  of  k  with  X  arbitrary  and  for  arbitrary  values  of  k  with 
X  =  1 .  These  solutions  were  used  to  check  the  numerical  procedure. 


h2~r'1/2v2(e) 


as  r-0  on  » 


(3.6UJ 


e^e^r.-e)  ,  e2  =  -e2(r,- e)  on  £  ,  (3.6i; 

where  pQ  is  a  positive  constant  and  \  a  constant  satisfying  (3.36). 
In  (3.59),  f  (e)  is  3  smooth  function  governed  by  (3.56),  (3.57).  For 
m=-1/2  the  implicit  equation  for  »p( e )  (3.18)  can  be  inverted  analyti¬ 
cally  as  in  the  anti-plane  shear  problem.  This  result  and  (3.14)  even¬ 
tually  yields 


v-,(e)  •  1?C°S  ’ si"  9  '  (3'62) 


0  (cose  +  P(e))£ 


where 

p(e)=  (4cos2e-3)1/2  (-J<e<f)  • 

Further,  (3.13),  (3.24),  (3.25)  and  (3.62)  provide 


(-f<e< 


7T  \ 


(3.63] 


(3.64; 


Alternative  representations  for  v^  are  obtained  from  (3.29),  (3.49) 
and  (3.51)  in  the  neighborhood  of  e  =  0: 
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v-, -fPo1/3t1  +j92  +  0<94)]  ' 

►  as  0  -0  .  (3.65) 

v2  * 

Finally,  in  (3.59)-(3.61 ) ,  ?+  and  £  are  elliptic  regions  and  U 

the  hyperbolic  region  given  in  (2.20).  The  domains  are  separated  by  the 
elastostatic  shocks  S+  and  S".  Referring  to  (2.18),  (2.19),  (3.47), 
(3.51)  and  (3.52)  we  find 

1/3 

ST:  e  =  ±0(r)  ,  9(r)  ~  )  (|-  f.  (0))r1/2  as  r-0  .  (3.66) 


4.  HIGHER  ORDER  CONSIDERATIONS  AND  RESULTS 

What  remains  to  complete  the  first  order  analysis  is  the  evaluation 
of  the  constant  X  and  hence  the  functions  7  (e)  on  [0,u].  In  this 
section  we  calculate  these  unknowns  numerically  through  higher  order 
considerations  and  the  subsequent  results  are  presented. 

To  begin,  we  replace  (3.66)  by  the  two  term  asymptotic  representa¬ 
tion 

S_:  e  =  ±©(r)  ,  o(r)  ~  {--)  (f-  7,  (0))r1//2  +  Br*  as  r  —  0  ,  (4.1 ) 

where  B  and  t>l/2  are  undetermined  constants.  The  domains  *,  s+ 
and  C  remain  as  defined  by  (2.20).  Similarly,  (3.60)  is  modified, 
such  that 

h.j  ~r3^2v-j  ( e )  +  r  ^(e) 

►as  r-0  on  V  .  (4.2) 

h2~r'1/2v2(e)  +  r  ^(e)^ 

In  (4.2)  n^  >3/2  and  n2>-l/2  are  constants  while  w^(e)  and  w2(e) 
are  smooth  functions  on  ( -tt/6 ,tt/6 ) .  Equations  (3.27),  (3.49)  and  (3.51) 
provide  the  Jacobian  to  leading  order.  Now  assume 

J  ~r"^p2^(e)  +  r  ^  (e)  as  r-*0  on  V  .  (4.3) 


The  constant  is  assumed  to  be  greater  than  -1  and  ( e )  has  the 
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same  smoothness  as  p ( e )  on  ( -ir/6 ,tt/6 ) ;  see  (3.13). 

Observe  that  (3.49)  and  (3.51)  imply 

R~r4u(e)  as  r-0  on  U  ,  (4.4) 

with  u(e)  a  continuously  differentiable  function  for  -ir/6<e<ir/6. 

An  alternative  asymptotic  representation  for  R  is  generated  from  the  defini¬ 
tion  (3.4)-(3.6),  together  with  (4.2)  and  (4.3).  On  equating  this  ex¬ 
pression  to  (4.4)  we  arrive  at 

v,+l  _oio  n«+ 1/2  ?!  a 

3r  p”*7  q-j- 2r  c  p“^(v2w2  ~  2"n2v2w2 ^  ~r aS  *  (4-5) 

on  U.  Suppose 

v1*n2-l/2<3  ;  (4.6) 

then  (4.5)  and  (4.6)  give 

3p4i/3q1  -  2(v2w2  -  ^-n2v2w2)  =  0  on  (-£,|-)  •  (4. 7) 

The  governing  differential  equation  (3.3)  with  a =2,  with  the  coordinates 
replaced  by  (4.2)  and  the  Jacobian  by  (4.3),  (4.6),  implies  that  w2>  q1 
also  satisfy 

p2/3w2  "  jP’1/3pwz  +  nz(n2  +  2)P2/3w2  -  2v2q-| 

+[v2+ (n2-|)v2]q1  =  0  ,  (4.8) 

for  -ir/6< e  <ir/6. 

It  is  convenient  at  this  point  to  assume  that  suitably  smooth  solutions 
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exist  to  (4.7)  and  (4.8)  on  (-ir/6,rr/6) .  For  small  |e|,  we  find  by 
making  use  of  (3.29),  (3.51)  and  (3.65)  in  (4.7)  and  (4.8),  that 


w2  =  p18+o(0  ) 


vK,/V°<e2> 


>  as  e-0 


(4.9) 


where  p-j  is  a  positive  constant  of  integration. 

Next,  we  consider  the  Jacobian  to  higher  order.  Through  (3.5), 
(4.2),  (4.3)  and  (4.6)  we  find 


n. —2  .  ^  ^  n«  o  .  . 

r  (?v2^1  +  ni^2wl )  ~r  (q^  - +  n^w^  as  r-0  ,  (4.10) 


on  V.  Take  n^  to  have  its  smallest  possible  value 

nj  =  n2  +  2  , 


(4. II)1 


whence 

v2w-|  +  2(n2  +  2)v2w^  =  Hq^Sv^Wg  +  2n2V|W2  on  (-£>£)  •  (4.12) 

Again  we  assume  the  existence  of  a  smooth  solution, while  for  small  |e| 
equations  (4.9)  and  (4.12)  give 

P"4/3P 

V  3^ynT+0(e2>  as  9  -°  •  <4-,3> 

Now  turn  to  the  solutions  of  the  differential  equations  (3.3)  on 
^Note,  if  ni <n2  +  2,  wi  is  unbounded  on  e  =  0. 
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the  elliptic  domain  s+.  We  replace  (3.59)  by 


e.j  ~  x 


r  sin  9  +  p 


(e) 


1  2 


1  3 


e9~x^pj  ?i  (0))  -  xr  cos  e  +  p„  3r2  f?(e) 


>  as  r-0  on  s 


(4.14) 


_  2 

+  poTr2g2(e) 


where  ga(e)  are  analytic  functions  on  [0,ir],  In  (3.3)  we  assume  the 
coordinates  conform  to  the  representation  (4.14)  and  obtain  the  follow¬ 
ing  differential  equations  for  g^: 


La(gl  ,92;2,A)  =  Ha(fl  »f2’X)  for  °<0<1T  »  (4.15) 


where  the  differential  operators  La  are  defined  in  (3.39)  and  (3.40). 
Furthermore, 


H1(f1,f2;x)  =  4x  ^[(J1T1 )‘+  3J1f1]  -  6X  T(sineJ1  +  2cos  eJ,  )J] 

4  _7  _  7 

+  x^[3R1F2  -  ^^2  + 2sin  0(2X  3J2  -  XR2)  +  2cos  0(2X  3J2-xR2)‘] 

_  2  _  1_ 

H2(fi  ,?2;x)  =  4x  3[(0^72)’+ 3J-|T2]  -  6x  3(2sin  0O1  -  cos  eJ^  )0-| 


(4.16) 


1  .  .  .1  .1  (4J7) 

-  x3[3R1F1  -  R]F1  +2cos  0(2X_1J2  -  x  3R2)  -  2sin  e(2x']  J2  -  x  3R2)’], 
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on  (0,ir)  with 

-4 .  .  .1 

J1  =  x[sine(|f1  +  x  Tf2)  +  cos  0(7,  -fx  3?2)]  , 

r  (4.18) 

J2*  1^1*2  "  ^1*2^  ’ 


_  4 

R1  =  A-1  (sin  e[|x  3(2  +  +  (2 +  M)?2] 

_  4 

+  cos  0 [ X  3(2  +  K)t1 -|(2  +  M)?2]]  ,  (4.19) 

8  _  8 

R2  a  -X"2{^(2  +  K  -  3x3cos  29 )f^  +  ^(2  +  M  +  3X  3cos  2e )f2 

8  _  8 
+  (2+  K+  3x3cos  2e)7^  +  (2  +M  -  3x  3cos  2e)?| 

8  _8 

+  9 ( x3 T|T^j  -x  3T2T2)sin  2e -3N[|(f1f2  + 7^2)00520 

-  (7i72 -|F172)sin  20]}  .  (4.20) 

In  (4.18)-(4.20)  the  constants  K,  M  and  N  are  supplied  by  (3.41). 

The  boundary  conditions  on  the  crack  face  0  =  tt  (3.8)  suggest 

.4  _  7 

3g1U)-2x  Tg2(ir)=|x  T(3f*(w)  -  f^(ir))  , 

_  4  _  7 

g2(")  -  2x  '^g1(ir)  =  -3X  ^7j  (Tr)72(ff) 


(4.21) 
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Associated  with  (4.15)  is  the  homogeneous  solution 


g^=G-][K+  (2+  K)cos  2e]  +  G2Msin2e  -  G3Ncos  2e 
g£  =  G4[M  -  (2  +  M)cos  20]  +  GgNcos  2e  +  G3Ksin  2e 


Y  (4.22) 


j 


on  [0,w],  where  G-j  -  G^  are  arbitrary  constants.  The  formulation  of 
the  boundary  value  problem  for  gQ  on  [0,ir]  is  completed  once  bound¬ 
ary  conditions  on  6  =  0  have  been  established. 

We  proceed  now  with  matching  the  higher  order  solutions  on  V  and 
C  across  the  elastostatic  shock  S+.  Recall  that  the  smoothness  of  the 
complete  solution  to  the  problem  (3.3)-(3.8)  requires  continuous  coordi¬ 
nates  and  Piola  tractions  across  S+. 

Equating  the  coordinates  on  using  (3.44),  (4.1),  (4.2),  (4.9), 
(4.11),  (4.13)  and  (4.14)  leads  to 


t- 


-Ar  ~p  Br 
o 


+  (f)  (r7i(0))Pir 


n2+? 


"a 


>  as 


n2+2 


r-  0 


(4.23) 


In  (4.23)  we  let 


n2  *  0  ,  t  =  1  ; 


(4.24) 


referring  to  (3.57)  suggests 
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1 

b--(4)j(§-Ti(0))pi  - 

2 
p  3 

P1  =  (j_  f°(Q))  [^X  V°^V  +  9l(0^ 

(4.24),  (4.25)  and  continuity  of  the  tractions  (3.45)  across  S+  reveals, 
after  considerable  algebra,  that 

_  4  7  _  8  1 

3^(0)  -2X  Jg2(0)  =  -ix-5  +  ix  3(X  3 -3)^(0)  -f  X3 ff(0) 

_  7  _  4  4 

+  |x_I^(0),  (^-71(0))?2(0)  +  []-x"3-2x34-71(0))]g1 

1  _8  _8  _i 

=  A3[]*x  3-^(9  +  25x  T)F-,(0)  +  (3  +  2X  3Jf^(0)]f2(0)  ,  J 

What  remains  is  to  solve  the  boundary  value  problem  (4.15)-(4.21 ) 
and  (4.26).  Recalling  (4.21),  (4.22)  and  (4.26)  we  note  that  the  homo¬ 
geneous  system  has  the  nontrivial  solution 

_  4 

g^(e)  =  2x  3 MXjSin  2e  , 

_  4  _  4  _  4 

g|(e)  =  £M(X  3  N  -  3M)  +  [2x  3  N+  (2  +  M)(3M -X  3N)]cos20jX1  , 

Vo r  X  obeying  (3.36),  l/2-71(0)^0. 
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where  is  an  arbitrary  constant.  Consequently,  in  order  for  a  solu¬ 
tion  to  the  nonhomogeneous  problem  to  exist  the  eigenfunction  (4.27) 
must  be  orthogonal  to  the  right  hand  sides  in  (4.15). 

The  particular  solution  to  (4.15)  is  found  in  terms  of  A  and  f 
through  the  method  of  variation  of  parameters.^  This  expression,  with 
(4.22),  facilitates  a  numerical  investigation  into  the  existence  of  a  so¬ 
lution  to  the  boundary  value  problem.  We  find  functions  g£(e),  satisfy¬ 
ing  the  nonhomogeneous  differential  equations  (4.15)-(4.20)  and  the 
boundary  conditions  (4.21),  (4.26),  if 


3 

A  =  A  =  2. 657 <  (2  +  t/J)4  (=  2.685)  .  (4.28)2 

This  result,  together  with  the  eigenfunction  (4.27),  provides  the  com¬ 
plete  solution: 


ga(e) =  g^(e) + g^(e)  on  [O.n] 


(4.29) 


The  numerical  solutions  obtained  for  f  and  g*5  with  a  =  a  are  shown 

a  a 

Fig. 8.  Some  pertinent  results  from  the  figure  are 


f1  (0)  =  -  .065,  f2(0)  *  0,  f1  (it)  =  0,  f2(ir)  =  -  .090 


gP(0)*-.020,  gP(0)  =  0,  gP(»)=  .105,  gg(*)i-.152 

Recall  that  the  functions  fa  on  [0,ir]  depend  uniquely  on  A  for 
Aj*A*  given  in  (3.58).  If  A  *  A*,  7a  comprises  of  an  eigenfunction 
and  a  particular  solution. 

2 

The  result  (4.28)  Is  confirmed  by  an  independent  numerical  computation 
in  which  (4. 1 5)- (4 . 20)  subject  to  (4.21),  (4.26)  were  treated  directly 
using  a  standard  differential  equation  solving  routine. 


-  (4.30) 
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>0  for  -J<9<J  , 

yv|(9) -iv|(e)4  (4.33) 

=  0  for  e  *  ±  g-  , 

and  thus  is  a  smooth  function  on  (- ir/6,ir/6).  Furthermore,  (4.12) 
and  (4.24)  provides 

l0l  ,  _i  V 

Wi(0)  =  V24(0)  v^vMyP  ^(«p)Vg(cp)  "  3vi »  (4.34) 


which  also  has  the  appropriate  continuity  on  (- tt/6,w/6) .  The  constant 
Pj  is  evaluated  through  (4.25),  (4.28)  and  (4.30)  to  reveal 


035 


(4.35) 
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which,  in  conjunction  with  (4.9),  (4.13)  and  (4.24),  gives 


W1  =WlP-3+°(e2),  =  .006 


>  as  e  -  0  . 


w2  =  W2pqJ  9+0(e3),  W2=  -  .035 


(4.36) 


Summary  of  the  final  results 

From  Eqs.  (4.2),  (4.11),  (4.14),  (4.24)  and  (4.28)  we  observe  that 
the  second  order  asymptotic  solution  to  the  local  crack  problem  (3.3)- 
(3.8)  is  as  follows: 


3 

—  —  _o 

yj  ~djr  sin  9  +  rtj  (e)  +7^(9) 


y2~c2  +  d2r cos  0  + ^^(e)  +  r2g2(e) 


>  as  r-0  on  £ 


(4.37) 


_  -2— 
y1  ~r  v1  (e )  +  r  w1  (e) 

1 

—  —~7—  — 

y2~r  v2(e)  +  w2(9) 


fas  r-0  on  V 


j 


(4.38) 


where 


(4.39) 


and  the  coordinates  on  5  are  obtained  through  symmetry.  The  constants 

c9,  d  are  found  to  be 
c  a 
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c2=  1.013,  d]  =  .722,  d2  = -  2.657  .  (4.40) 

Recall  that  ga(e) =  g^(e) + g^(e)  on  [0,n]  in  (4.37).  An  analytical 
expression  for  g^e)  is  given  in  (4.27)  while  g^e),  together  with 

a  a 

fa(e),  are  described  in  Fig. 8.  In  (4.38) 

1 

V0)  =  po  vl  (Q) »  ^(0)  =  Pq1v2^6^  ’ 

2  _  2 
if,  (e)  =  p3  w1  (e),  w2(e)  =  po3w2(6)  , 

where  v^(e)  and  v2(e)  are  found  in  (3.62)  and  (3.64),  respectively. 

We  refer  to  Eqs.  (4.31 )-(4.36)  to  obtain  expressions  for  wa(e).  In  ad¬ 
dition,  we  note  that  the  domains  and  U  (Z.20)  have  the  elastostatic 
shock  S+  (2.18),  (2.19)  as  a  common  boundary.  In  view  of  (4.1),  (4.24), 
(4.25),  (4.28),  (4.35)  and  (4.39)  the  asymptotic  representation  of  the 
shock  reduces  to 

1 

S+:  9~c2r^  +  c jF  as  r-*0  .  (4.42) 

The  constant  c2  is  given  in  (4.40)  and 

c3=.036  .  (4.43) 

In  the  analysis  sumnarized  above  the  remaining  unknown  constants  are 
P0  (>0)  and  *i  in  gQ  (see  (4.27),  (4.29)). 

The  results  (4. 37)- (4. 43)  represent  a  locally  one-to-one  map  of  the 
neighborhood  of  the  crack-tip.  Calculating  the  deformation  image  of  the  crack 
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face  e  =  it  in  terms  of  the  nondimensional  ized  spatial  coordinates  leads  to 

1 

y2~Vvf  (c2=  1.013,  =  8.24)  as  y^O*  .  (4.44) 

The  crack-tip  (r=0)  deforms  onto  the  line  -c2sy2sc2*  yl  =  0  and  the 

+  *4* 

shock  S  (4.42)  maps  onto  S  which  is  found  to  satisfy 

2 

y2  ~c2  ■  c5^1  (c5  *  3.48)  as  y1-0+  .  (4.45) 

The  curves  (4.44),  (4.45)  are  sketched  in  Figure  9  where  V  and 

*+  + 

£  correspond  to  the  respective  images  of  X  and  s  .  It  is  apparent 

from  the  figure  that,  in  the  vicinity  of  r=0,  the  traction  free  faces 
(e  =  +ir)  open  to  almost  flat  surfaces  oriented  in  the  direction  of  the 
applied  load  at  infinity.  The  fact  that  the  curves  extend  slightly  into 
the  half-plane  y-j  > 0  cannot  be  explained  by  this  local  analysis. 

We  now  list  the  Cauchy  stresses  as  functions  of  the  material  coordi¬ 
nates  (r,e)  in  the  upper  half-plane  computed  from  (1.30),  (1.33)  and 
(4.37),  (4.38): 


T11 

on 

C+  , 

■—  —  ■* 

V 

o(r4) 

on 

V  , 

"(1  -  a"8/3)  =  .926 

on 

C+  , 

t22 

-  -i 

V 

1 

on 

► 

X  , 

T1 2 

V 

on 

+ 

-  rs*  ■* 

U 

o(r2) 

on 

x  , 

where 


1 

o  ^  *a  _  i  ^  o  _ 

(e)  =  3a  (cos  ef^  +  |-si:  af  j )  +A  (sin  ef2  -  -j-cos  ef2) 

7 

E12(e)  =  A"  (sin  ef-j -■g-cos  ef-j) +A  (cos  ef2 +  £Sin  ef2) 


^  (4.47) 


j 


and  recall  A=  2.657.  In  the  hyperbolic  region  V  the  leading  terms  in 
the  coordinates  do  not  contribute  directly  to  the  stresses  and 

t-j 2 .  The  functions  E^(e)  and  E12(e)  on  tO»*]  are  plotted  in  Fig. 10. 
The  stresses  (4.46)  indicate  that  the  material  in  the  neighborhood  of  the 
crack-tip  is  subjected  to  deformation  that  approximates  uniaxial  tension 
parallel  to  the  y2-axis. 
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APPENDIX  A 


In  this  Appendix  solutions  to  the  anti-plane  shear  crack  problem 
are  found  which  include  two  elastostatic  shocks  that  are  asymptotically 
tangent  to  rays  issuing  from  the  crack-tips. 

To  this  end,  let  S+  and  S”  be  curves  defined  in  (2.18),  (2.19) 
with  s  =  0  and  A=e*,  so  that 

ST:  e  =  ±e(r)  ,  0<rsrQ  ,  (A.l) 


where 


0(r)~9*  as  r-0  .  (A. 2) 

The  regions  and  described  in  Figure  11  are  defined  locally 

as  follows 


*"■)  s  C(r,0)|  -  9(r)<e<0(r)  ,  0<rsrQ3  , 
^•£(r.e)|e(r)<e<,  ,  0<rsrQ3  , 
Ur»e)l  -  ir<e<-©(r)  ,  0<rs rQ)  . 


>  (A. 3) 


Keeping  in  mind  the  synmetry  imposed  by  the  deformation  at  infinity  (2.3) 
we  assume  the  displacements  admit  the  representation 


4 
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m, 

u~r  v^(e)  ,  v-|  (e)  = -v-|  (-e)  as  r~0  on  ^  ^ 


m. 


2  + 
u~r  V2(e)  as  r-0  on  » 


u~-r  (-0 )  as  r- 0  on 


>  (A. 4) 


J 


where  the  exponents  ma  satisfy 


0<m  <1  .  (A. 5) 

a 

Furthermore,  v^(e)  and  V2(e)  are  smooth  functions  on  the  overlapping 
intervals  (-9^,9-j)  and  [ir,92)»  respectively  -  the  angles  e1  and  92 
(yet  to  be  determined)  obey  0*£0j  a  it  and  O<e2*0*.  The  boundary  con¬ 
ditions  (2.4)  imply 

v2(w)®0  ,  (A. 6) 

We  repeat  the  analysis  resulting  in  Eqs.  (2.10)-(2.15)  for  the  pres¬ 
ent  case  and  find  the  general  solutions  for  va(e): 

m  v  (0 )  =  V  cos  ( © )  |  v  +  cos  2if/  ( 0 ) {  ma^  ,  (A. 7) 

a  a  a  a  1  a  a  1 

where  \i>  (e)  is  the  solution  of 
a 

(i+  cos  2tp  )<ji  +  (v  +  cos  2ip  )  =  0  .  (A. 8) 

'3  Tar  a  a  a 

The  constant  V  is  assumed  positive  while 
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Vj(3-2ma>  •  (A-9> 

We  note  also  that  v  -t-  cos  2^  /0  in  order  that  the  displacements  be 

a  a 

bounded. 

In  what  follows,  we  consider  Ose^ir  and  refer  to  (A. 4)  for  the 
displacement  in  the  lower  half-plane.  As  in  (2.16),  let 

(0)  -  w/2  ,  (A. 10) 

so  that  v-| (0)  =  0.  The  boundary  condition  (A. 6)  is  satisfied,  without 
loss  of  generality,  by  stipulating 

^(ir)  =  0  .  (A. 11 ) 

The  solutions  ^(e)  and  ^(e)  governed  by  (A.8)-(A.ll)  are  sketched 
in  Figure  12.  The  function  ^(e)  decreases  monotonically  from  a  value 
of  tt/2  at  e  =  0  (as  specified  by  (A. 10))  and  approaches  the  asymptote 
ij7=  1/2  cos'1  (-v-j )  in  the  limit  e-»oo.  ^(e)  is  a  decreasing  single 
valued  expression  on  the  interval  [eo,ir]  where 

( 1  -  nip)  -I  i 

eQ  =  *  -  V - —  tanh"  (-*tan  *  )  , 

'  mgv  v 

(A. 12) 

1 

1  ,  Jrnc-1/  h 

,  ^-^cos  (-j)  . 

Further,  we  note  that  wh,,1e 

On  substituting  <1^(9)  into  (A. 7)  we  find  v-j(e)  and  V2(b)  are 
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smooth  functions  defined  on  [0,*)  and  (0o»tt]  respectively  -  the  in¬ 
terval  (0O.*)  is  common  to  both  domains  in  accordance  with  the  original 

assumption  (where  8,=*  and  0,=  0j.  Also  V  >0  implies  v  have 

I  do  a  a 

the  correct  sign,  namely 


v-|(0)>0  for  0<9<tt,  V2(0)>O  for  tt  ss:  e  <  eQ  . 
Accordingly,  the  shock  angle  0*  must  satisfy 


(A. 13) 


0  <  0 "<  7r 
0 


(A. 14) 


An  asymptotic  equilibrium  solution  displaying  adequate  smoothness  is 
now  generated  to  this  anti -plane  shear  crack  problem  by  matching  the  dis¬ 
placements  and  associated  tractions  across  the  shock  S+. 

Continuity  of  displacements  (A. 4)  requires 


m-i  m«  . 

r  Vj(e)~r  v2(e)  as  r-0  on  s 


(A. 15) 


On  referring  to  (A.l),  (A. 2)  and  (A. 13)  we  conclude  that  asymptotic  equal¬ 
ity  (A. 15)  holds  if  and  only  if 


i*!  =m2  =  m,  v1  (0*)  =  v2(0*) 


(A. 16) 


This  result,  together  with  (A. 7)  provides  one  matching  condition  across 


the  shock: 


v  +  cos  2^(0*)  "  ^O-m)  cos<|>2(0*) 
V2  v  +  cos  2^(9*)  cos  1^(0*) 


(A. 17) 


-60- 


where 

w- v  =1(3-  2m)  (l<v<1)  .  (A. 18) 

a  J  J 

Further,  (A. 14)  Indicates  that 

^<i»1(e*)<J  ,  0<^2(e*)<^o  (4>0<i)  .  (A. 19) 

The  remaining  condition  to  hold  on  S+  pertains  to  the  matching 
of  tractions.  Equations  (1.18),  (1.24),  (1.25),  (A. 4)  and  (A. 7)  imply 

-1/2(1  -m)  sin2tp,  (©*) 

=  - - -  .  (A. 20) 

sin  i|>2(6*) 

In  view  of  (A. 17),  (A. 20)  is  replaced  by 

p  p 

sin  (e*)cos  i|)-|  (e*)  =  sin  i|>2(e*)cos  <|i2(e*)  ,  (A. 21) 

which,  in  turn,  can  be  written  in  the  form 

♦1(0*)-^2(0*))  for  0<<|>2(e*)<^0  *  (A. 22) 

Here,  3(*)  is  a  monotonically  decreasing  function  on  [0,^]  with  the 
endpoints  3(0)  =  ir/2  and  ^q)  =  i>0- 

What  remains  is  to  prove  the  existence  of  a  solution  to  (A. 22)  for 
an  admissible  value  of  0.  The  shock  condition  (A. 17)  can  be  made  to 
hold  merely  through  an  appropriate  choice  of  the  ratio  V-j/Vg.  From  the 
sketch  of  3(ij>2(e))  In  Fig. 12  It  Is  apparent  that  we  can  find  an  angle 
0*  (depending  uniquely  on  m)  for  which  the  curves  4^(0)  and  3(i»2(e)) 


v  +  cos  2^1(0*) 


v  +  cos  2g»2( 0*) 
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intersect  and  thus  satisfy  (A. 22). 

Consequently,  there  exists  a  system  of  solutions  to  the  crack  prob¬ 
lem  (2.2)-(2.4)  in  which  the  displacement  field  near  the  crack-tip  is  of 
the  form 


u~rmv(e)  as  r-*0  (-usesit) 


(A. 23) 


where 


V  (  0 ) 


^v-|(9)  for  -e*<e<6*  , 
v2(e)  for  e*<0^w  . 


(A. 24) 


In  (A. 23)  the  exponent  0<m<l  remains  undetermined  from  the  leading 
order  analysis.  The  function  v(0 )  on  [-ir.ir]  (A. 24)  possesses  a  dis¬ 
continuous  first  derivative  across  the  elastostatic  shock  S+  (A.l), 

(A. 2).  The  shock  angle  0*  and  v(e)  depend  on  the  unknown  exponent 
m  while  v(e)  is  also  arbitrary  with  respect  to  a  multiplicative  con¬ 
stant.  An  appeal  to  symmetry  reveals  the  displacement  on  Further¬ 

more,  the  displacement  equation  of  equilibrium  is  everywhere  hyperbolic 
as  r-0  for  -irseiir. 

We  now  consider  the  stresses  associated  with  the  deformation  (A. 23). 
The  nontrivial  components  of  the  actual  stress  field  are  obtained,  in 
terms  of  the  polar  coordinates  (r,0,z)  =  (r,0,Xj),  from  (1.18)  and  (1.19): 


t  =  2W'(  I  )  —  — 
T©z  "^Vr  30 


tzz  =  2W/(I1)|vu|2  , 


(A. 25) 


where  I^  =  3+|vu|  .  Referring  to  (1.24),  (A. 7),  (A. 8)  and  (A. 23)  yields 


i(l-m)  i-(l-m)  -  i 

Trz~yr  cos*(e)p  (e),  T0z~ur  simj;(0)p  (e)  , 


in-1  i 

x  ~yr  p  ( e )  as  r-0  (-irsasiO  , 


(A. 26) 


p(e)  *  V|v  +  cos  2ip(e )  | 


-JO-n) 


(A. 27) 


In  (A. 26)  and  (A. 27) 


*■,(6)  for  -0*<e<e* 

i|»(e)  =-.  i|»2(e)  for  e*<e*ir  ,  (A. 28) 

ir— 4^2  ( — ® )  for  -its0<-0* 

Finally,  V=V-j  on  (-e*,e*)  and  V  =  V2  on  (0*,tt],  [-it , -e*)  such 
that  the  ratio  /V2  is  given  in  (A. 17) 

The  asymptotic  solution  presented  in  this  Appendix  was  not  found  in 
[4].  In  that  study  the  applied  shear  at  infinity  y  was  assumed  small 
compared  with  one  and  the  solutions  to  the  appropriate  displacement  equa¬ 
tion  of  equilibrium  in  the  vicinity  of  the  crack-tip  where  chosen  so  as 
to  facilitate  matching  onto  the  linear  elasticity  solution  valid  else¬ 
where.  We  note  that  the  hodograph  transformation  used  in  [4]  to  solve 
the  small-scale  problem  may  not  yield  all  possible  solutions.  Further¬ 
more,  for  more  severe  deformations  at  infinity  in  which  y  is  not  neces¬ 
sarily  small,  the  results  in  [4]  are  no  longer  appropriate.  In  such 
cases,  displacements  of  the  form  (A. 23)  may  occur  near  the  crack-tips. 
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APPENDIX  B 

Here,  we  establish  that  no  smooth  solutions  of  the  form  (3.9)  exist 
for  the  problem  (3.3)-(3.8).  Accordingly,  we  let 

m 

y  ~r  “v  (e)  as  r-0  (no  sum)  ,  (B.l) 

a  a 

uniformly  for  -ir^esrr,  where  are  constants  such  that 

Oimin[m^  ,nu]<l  .  (B.2) 

Further,  we  stipulate 

va*0  ,  (B.3) 

are  smooth  functions  on  [-ir,ir].  The  symmetry  (3.7)  associated  with  this 
problem  implies 


v-j (0)  =  v1  (-0)  ,  v2(e)  = -v2(-e)  .  (B.4) 

The  inequality  (B.2)  assures  singular  deformation  gradients  and  bounded 
displacements  at  r=0. 

Suppose  m^=m2s0  in  (B.l),  whereupon  the  coordinate  equations  of 
equilibrium  (3.3)  are  satisfied,  to  leading  order,  if 

p^(e)  =  ( 0 )  +  Vg ( 0 )  =  0  on  [-it, u]  ,  (B.5a) 


or 


l 
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J~Qr"^^3p'(e)  ,^(0)^0  on  [-tt,tt]  as  r-0  ,  (B.5b) 

where  Q  is  a  positive  constant.  Equations  (B.3),  (B.4)  are  not  con¬ 
sistent  with  (B.5a)  on  [-tt.tt].  Now  assume  (B.5b)  holds  and  note  that  equi¬ 
librium  of  an  arbitrary  region  js=[(r,e)  |O<r<ro,0^<e<02»6i<02  e  (-tt.it)} 
requires,  in  the  limit  as  ro-0,  that 

p^(©)  =  Q  ^  on  (-tt.it)  .  (B.6) 

According  to  the  boundary  conditions  (3.8), 


P^(9)  =  0 


on  0  =  -it, it 


(B.7) 


The  results  (B.6)  and  (B.7)  are  incompatible  with  the  assumed  smoothness 
of  va  and  thus  eliminating  the  case  m1=m2  =  0. 

Now  let 

m^=m2  =  m  ,  0<m<l  ,  (B.8) 


in  (B.l).  Replacing  the  coordinates  in  (3.3)  by  (B.l)  and  retaining 
dominant  powers  of  r,  provides  two  coupled  nonlinear  differential  equa¬ 
tions.  Assume  the  equations  hold  on  (-tt,tt)  and  let 

v-j  (0)  =  p(o)cos  9(0),  v2(e)  =  p(e)sin 9(0)  on  [-tt.it]  .  (B.9) 


If 


n(e)  s  p(fl)/p(e),  x(e)  =  <p(&)  on  [-tt.tt]  ,  (B.10) 


the  differential  equations  reduce  to 


r 
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(3m2  +  n2)x  -  [nn  -  m(m  +  4)n  -  n^]x  +  nx^  =  0  »  (B.ll) 

2mnx  -  [mn  +  m3  -  (4  -  m)n2]x  +  (4  -  3m)x3  =  0  ,  (B.12) 

on  (-7f,ir).  The  boundary  conditions  on  the  crack  faces  (3.8)  imply 

v2(e)  +  v2(e)  =  0  or  v2(e)  +  v|(e)  =  0  on  6=-*,*  .  (B.13) 

In  view  of  (B.9)  and  (B.10),  (B.13)  is  equivalent  to 

p(e)  =  0  or  x(e)  =  0,  p(e)  =  0  on  e  =  -ir,ir  .  (B.14) 

Solutions  of  (B.ll)  and  (B.12)  subject  to  (B.14)  can  be  examined  in 
the  phase  plane.  For  convenience,  let 

ft(0)  =  n^(8)»  X(e)  =  x^(e)  on  ,  (B.15) 

and  derive,  from  (B.ll)  and  (B.12),  the  trajectory  equation 

E'Tffef1  <x>0-n>0>  • 

with 

F(a,X)«  2(1  -m)X- [m2(2  +  m)  -  (Z-m)fl]  ,  (B.  17) 

G(fi,X)  =  [3m2(4- 3m)+ (4- 5m)fljX  + [(4  -  3m)fi- 3m3](m2  +  fl)  .  (B.  18) 

The  trajectories  X=X(n)  are  sketched  in  Figs.  13  and  14.  From 

p 

the  figures  we  note  that  if  X(e*) =  x  (e*) =  0  (fi(e*)?*fl0)  for  some 
-irse*sir  then,  X(e)=0  on  [-*,*].  Further,  a  solution  to  (B.ll) 
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and  (B.12)  such  that  n  and  X  lie  on  the  separatrix  can  only  attain 

the  values  at  the  equilibrium  point,  q  =  q0  and  X=0,  in  the  limit 

e-*-oo  or  +co.  Consequently,  the  boundary  condition  (B.14),  in  which 

x(e)  =  0  on  9  =  -ir ,ir ,  can  only  be  satisfied  if  x(e)  =  9(0)  =  0  on  [-ir,n]. 

From  (B.9)  we  conclude  v  (e)  =  c  p(e)  on  [~iT,ir]  while  (B.4)  implies 

a  a 

the  constants  c  =0. 

a 

Assume  now  the  first  condition  in  (B.14)  holds  and  consider,  in  par¬ 
ticular,  the  crack  face  0 = w.  From  (B.10)-(B.12)  and  (B.16)-(B.18)  we 
obtain  an  implicit  representation  for  p(e): 

p(e)  =  eJ(0)  (B.19) 


on  C-ir.Tr],  with 


J(e) 


X(e) 

X(tt) 


m(3m2  -  n(X)) 
4XF(g(X),X) 


dX 


(B.20) 


In  (B.20),  the  appropriate  trajectory  n  =  n(X)  (X>0)  is  chosen  with 
regard  to  the  boundary  condition  p(n)  =  0,  which  dictates 

<4(8)  -  -oo  as  9-*ir  .  (B.21) 


We  investigate  the  possibility  of  satisfying  (B.21).  Suppose  T 
is  a  curve  n=fl(X)  in  the  phase  plane  that  intersects  the  straight 
line  F(fi,X)  =  0  when  9  =  ir  (see  Fig. 13).  Thus,  the  integral  (B.20) 
is  singular  at  the  boundary.  Further  analysis  indicates  the  integrand 
is  0[(X-XU)r1/2]  as  X-X(ir)  on  T  which  implies  the  condition 
(B.21)  is  not  attained.  Now  let  X(it)  =  oo.  It  can  be  shown  from 
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(B.16)-(B.18)  that 


4- 5m) 

TT^rj 


Q  =  0\X 
n  =  0(log  X) 
n  =  o(1) 


(0<m< 4/5)  ^ 
(m-4/5) 
(4/5<m<  1 )  j 


>as  X-oo 


(B.22) 


Using  (B.22)  in  (B.20)  we  again  find  the  integral  remains  bounded  for 
0<m<l.  It  is  apparent  then,  that  suitably  smooth  solutions  of  (B.ll) 
and  (B.12)  cannot  be  made  to  satisfy  (B.14). 

What  remains  is  to  account  for  the  case  Accordingly,  let 

m2  =  m<mi,  0<m<1  .  (B.23) 


Define  the  function  q(e)  on  [-it, it]  through 


J~rvq(e)  as  r-0  ,  (B.24) 

where  (B.l),  (B.23)  and  (3.5)  indicate  that  the  constant  v>2(m-l). 
Further,  we  assume  q  is  continuously  differentiable  on  (-tt.tt)  and  is 
continuous  up  to  the  boundaries  9  =  -tt,tt. 

Substituting  (B.l),  (B.23)  and  (B.24)  into  the  coordinate equil ibrium. 
equations  (3.3)-(3.6)  and  considering  only  the  leading  order  terms  we  ob¬ 
tain  two  nonlinear  differential  equations  for  q  and  Vg.  Eliminating 
q  from  one  equation  yields 


1  Note  that  0<  (4-5m)/4(l-m)  <  1  for  0<m<4/5. 


i  .jr-'al.tafj— A.. 
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mv2(4m2v2  -  p2)(v2  +  iri2v2)  -  2(m- l)p2(2m2v2  + p2)  +  3vp4  =  0  ,  (B.25) 


(2m2v2  +  p2)q-^(2vmv2  +  p2)q=  0  , 


(B.26) 


for  -Tr<e<ir,  where 


p2(e) » m2v2(e) + v2(e)  , 


(B.27) 


on  [-Tr,?r] .  Equations  (B.l),  (B.23),  (B.24)  infer  the  boundary  conditions 


r3vq3(e)~r2(m‘1)m2v2(0)  as  r-0 
v2(e)v2(0) =  0 


on  0  =  -it  ,ir 


(B.28) 


The  technique  adopted  to  solve  the  differential  equation  (2.12)  is 
applied  to  (B.25).  We  obtain 


e,  e« 

mv2(e)  =  V2|cos  i^(e)|  |e0  +  cos  2ip(e)|  cos^(e)  ,  (B.29) 


e3 

q(0)-Q|cos*(e)|  |  e  +  COS  2i|>(e)  j 


(B.30) 


on  [-ir.ir],  where  V2  and  Q  are  positive  constants  and  t/*(e)  satisfies 


(^■+  COS  +  (eQ  +  cos  2ip )  =  0  , 


on  (-it, it).  The  constants  eQ  -  e3  are  given  by 


(B.31 ) 
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e0  =  2-f(m-v)  ,  £]  =  2(1 -eo)  _1 

m(1  -  2eQ)  2(v  +  1 )  (1  -  eQ)  -  m 

e2=  4(1  -  e0)  *  e3= 


Me0^l)  , 


(B.32) 


for  0<m<l.  We  note  that  if  0<eo<]  any  solution  ip(d)  of  (B.31) 
on  the  finite  interval  [-n-,7r]  must  necessarily  satisfy 


.  ♦1  =^cos_1(-eo)  (0<^<|)  .  (B.33) 

This  condition  is  violated  only  in  the  limit  e--oo  or  +oo. 

Equation  (B.29)  accommodates  the  smoothness  requirement  and  the 
parity  condition  (B.4)  if  we  choose 

i|/(0)  =  ir/2  and  e-j  =  0  or  e-|  *  1  .  (B.34) 

Taking  =  0  in  (B.32)  gives  e0=l-m/2  and  v  =  2/3(m-l).  The  bound¬ 
ary  conditions  (B.28),  together  with  (B.29),  (B.30)  and  (B.32),  then 
imply  without  loss  of  generality,  that 

ij)(ir)  =  0  ,  4>(-ir)  =  ir  and  0^  =  vjj  .  (B.35) 

Note  that  for  0<m<l,  1  /2 <  eQ <  1  and  assume  the  existence  of  a  smooth 

solution  to  the  boundary  value  problem  (B.31),  (B.34)  and  (B.35).  Ac¬ 
cordingly,  for  some  value  of  e  on  (-ir,ir)  <f»  must  attain  the  value 
4»i  ®  1/2  cos  ^(-eQ)  and  thus  contradicting  the  result  (B.33).  We  con¬ 
clude  then,  that  no  smooth  solutions  exist  for  <|>  and  hence  v^  on 
[-ir, it].  A  similar  conclusion  is  reached  for  the  case  e^l. 
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FIGURE  I.  RESPONSE  CURVE  IN  SIMPLE  SHEAR 
FOR  THE  SPECIAL  INCOMPRESSIBLE 
MATERIAL  CHARACTERIZED  BY  ( 1 : 23 ) 


A 


3 


4 


FIGURE  2.  BLATZ-KO  MATERIAL.  PLANE  STRAIN 
UNIAXIAL  STRESS:  r)faO .  NOMINAL  AND 

ACTUAL  STRESS  VS.  AXIAL  STRETCH 


FIGURE  3.  GEOMETRY  OF  THE  GLOBAL  CRACK 
PROBLEM 


I 
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FIGURE  4,  SKETCH  OF  THE  PHASE  ANGLE 
*(0)  VS.  0  FOR  0<  m<  I 


FIGURE  5.  UNOEFORMED  CONFIGURATION 

near  the  right  crack-tip 


5- 


FIGURE  6.  (m.k)-DIAGRAM  FOR  THE  ANTI 
PLANE  SHEAR  CRACK  PROBLEM 
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FIGURE  7.  (m,k)-DIAGRAM  FOR  THE  PLANE 
STRAIN  CRACK  PROBLEM 
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7„(9)  AN0gJ<9)  :  X=A 


0 


FIGURE  9.  DEFORMATION 
STRAIN  CRACI 


FIGURE  II.  LOCAL  GEOMETRY  NEAR  THE  CRACK-TIP 
FOR  THE  ANTI-PLANE  SHEAR  PROBLEM 
TREATED  IN  APPENDIX  A 
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FIGURE  12.  SKETCH  OF  THE  PHASE  ANGLE 

*  ( e )  VS.  9  REFERRED  TO  IN 
a 

APPENDIX  A 


FIGURE  13.  SKETCH  OF  THE  TRAJECTORIES  IN 
THE  PHASE  PLANE  FOR  EQUATION 
(B.  16):  0  <  m  <  4/5 


G(il,X  )s0 


/  _  3m  (4-3 m) 
2  (5m-4) 


FIGURE  14.  SKETCH  OF  THE  TRAJECTORIES  IN  THE 
THE  PHASE  PLANE  FOR  EQUATION  (B.I6) 
4/5  <  m  <  I 
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Cracks,  Blatz-Ko  material,  finite  deformations,  elastostatlc  shocks 


?n  TCny-^nTT^r^RW  WflHt?  deformations  of  an  elastic 

solid  Is  used  to  study  the  elastostatlc  field  near  the  tip  of  a  crack.  The 
special  elastic  materials  considered  are  such  that  the  differential  equations 
governing  the  equilibrium  fields  may  lose  elllptlclty  In  the  presence  of  suf¬ 
ficiently  severe  strains. 

The  first  problem  considered  Involves  finite  anti-plane  shear  (Mode  III)  defor¬ 
mations  of  a  cracked  Incompressible  solid.  The  analysis  Is  based  on  a  direct 
asymptotic  method.  In  contrast  to  earlier  approaches  which  have  depended  on 
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hodograph  procedures. 

The  second  problem  treated  is  that  of  plane  strain  of  a  compressible  solid 
containing  a  crack  under  tensile  (Mode  I)  loading  conditions.  The  material  Is 
characterized  by  the  so-called  Blatz-Ko  elastic  potential.  Again,  the  analysis 
Involves  only  direct  local  considerations. 

For  both  Mode  III  and  Mode  I  problems,  the  loss  of  equilibrium  ellipticlty 
results  In  the  appearance  of  curves  ("elastostatic  shocks")  issuing  from  the 
crack-tip  across  which  displacement  gradients  and  stresses  are  dlscontlnous. 
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